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ABSTRACT 


The  purpose  of  this  study  was  to  investigate  the  interaction  of  plane 
elastic  waves  with  a  thin,  hollow,  cylindrical  shell  embedded  in  an  elastic 
medium. 

The  cylindrical  shell  is  considered  to  be  elastic,  isotropic,  homo¬ 
geneous,  and  of  infinite  length.  It  is  surrounded  by  an  elastic,  isotropic, 
and  homogeneous  medium  whose  motions  conform  to  the  ordinary  theory  of  elas¬ 
ticity.  A  plane  stress  wave,  either  dilatational  or  shear,  with  a  step  varia¬ 
tion  in  time,  whose  wave  front  travels  in  a  direction  perpendicular  to  the 
cylinder  axis,  envelops  the  shell.  Later,  a  Duhamel  integral  is  used  to  study 
other  wave  shapes  for  the  incident  stress. 

The  response  of  the  shell  is  studied  by  expressing  the  two  components 
of  displacement,  radial  and  tangential,  in  terms  of  Fourier  series,  each  term 
of  which  is  called  a  mode.  The  equations  of  motion  of  the  shell  in  vacuo 
are  derived  from  expressions  giving  the  strain  and  kinetic  energies  due  to 
generalized  external  forces.  Forces  on  the  shell  result  from  the  stresses  in 
the  medium  at  the  boundary.  Stresses  in  the  medium  are  taken  to  be  the  sum 
of  the  stresses  due  to  the  incoming  stress  wave  expressed  in  terms  of  Fourier 
series  whose  coefficients  are  known,  and  those  due  to  the  reflected  and  dif¬ 
fracted  effects  expressed  in  terms  of  a  pair  of  displacement  potentials  rep¬ 
resenting  waves  diverging  from  the  axis  of  the  shell. 

The  equations  to  be  solved  consist  of  two  pairs  of  coupled  integro- 
differential  equations  in  the  generalized  coordinates  of  the  shell  and  the 
displacement  potentials.  By  use  of  a  digital  computer  they  are  solved  mode- 
wise  by  a  step-by-step  iterative  integration  technique  known  as  the  Newmark 
Beta  Method,  with  which  values  for  the  potential  functions,  and  the  acceler- 


ations,  velocities,  and  displacements  of  the  shell  are  determined*  Stresses 
in  the  shell  are  found  from  the  displacements,  and  the  values  of  the  potential 
functions  permit  determination  of  stresses  for  any  point  in  the  medium* 

Although  the  equations  are  written  to  include  an  infinite  lumber 
of  modes,  only  the  first  three  modes  are  considered  in  detail*  The  computed 
solution  is  compared  with  values  obtained  from  a  series  expansion  of  the 
equations,  which  is  valid  for  short  times,  and  with  the  static  solution  based 
on  the  theory  of  elasticity  to  which  the  general  solution  should  approach 
asymptotically.  In  addition,  the  results  of  two  particular  problems  are 
compared  with  results  given  in  another  study* 

Numerical  solutions  are  obtained  to  determine  the  effect  of  the 
several  parameters  which  describe  the  relative  physical  properties  of  the 
shell  and  medium.  Results  are  presented  in  tabular  and  graphical  form* 
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CHAPTER  I 
INTRODUCTION 


1.1  General  Remarks 

The  problem  of  designing  underground  protective  structures  to 
resist  the  effects  of  nuclear  weapons  has  become  increasingly  important  in 
recent  years  with  the  development  of  modern  weapons  whose  destructive  capacity 
is  overwhelming.  Engineers  in  this  field  are  hampered  to  a  great  extent  by  a 
lack  of  theoretical  information  on  how  structures  in  media  such  as  soil  or 
rock  behave  when  subjected  to  dynamic  loads.  Even  for  static  loads  alone, 
much  of  the  design  practice  today  is  of  a  semi -empirical  nature. 

When  a  nuclear  explosion  occurs,  stress  waves  are  transmitted 
through  the  air  and  ground.  How  are  they  transmitted  and  how  are  they  modified 
by  the  presence  of  a  structure  embedded  in  the  medium?  How  does  the  structure 
respond? 

The  purpose  of  this  report  is  to  study  one  aspect  of  the  problem, 
the  interaction  between  the  medium  and  structure. 

1.2  Statement  of  Problem 

The  problem  considered  here  consists  of  analyzing  the  elastic 
response  of  a  hollow  cylindrical  shell  (tunnel  lining)  embedded  in  an  elastic 
medium  when  subjected  to  an  incident  plane  stress  wave  traveling  in  a  direction 
perpendicular  to  the  axis  of  the  shell.  Some  questions  with  which  this  problem 
may  be  associated  are:  Do  tunnel  linings  in  contact  with  rock  afford  a  measure 
of  protection  significantly  higher  than  an  unlined  opening?  What  magnitude 
and  time  variation  of  displacement,  velocity,  and  acceleration  would  equipment 
mounted  within  such  a  structure  be  subjected  to?  How  are  stress  waves  within 
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the  medium  modified  in  the  vicinity  of  the  shell?  This  study  was  conducted 
in  an  attempt  to  find  some  qualitative  and  quantitative  answers  to  these 
questions  within  the  limitations  Imposed  by  the  assumptions  noted  below. 


1.3  Basic  Assumptions 

The  cylindrical  shell  is  considered  to  be  of  infinite  length,  and 
is  embedded  in  an  elastic  medium  of  infinite  extent  in  all  directions.  A 
plane  stress  wave  whose  front  travels  in  a  direction  perpendicular  to  the 
cylinder  axis  envelops  the  shell.  Strains  parallel  to  the  axis  in  both  the 
medium  and  shell  are  assumed  to  vanish;  thus,  since  each  cross  section  of  the 
shell  is  exactly  similar  to  every  other,  the  problem  is  reduced  to  one  of 
plane  strain. 

In  the  mathematical  development  of  the  problem  certain  basic 
assumptions  were  made,  the  most  important  of  which  are  given  here  with  a  few 
explanatory  remarks: 

(1)  The  medium  is  considered  to  be  homogeneous,  isotropic,  and 
linearly  elastic.  This  implies  that  the  ordinary  theory  of  stress  wave 
propagation  applies.  In  view  of  the  non -homogeneous,  non-isotropic,  and 
non-elastic  characteristics  of  most  materials  encountered  in  nature,  this 
is  a  severe  limitation;  however,  current  theories  of  stress  propagation 
through  such  media  have  not  advanced  to  the  stage  where  this  limitation  can 
be  readily  overcome.  In  the  case  of  some  rocks,  though,  this  assumption 
may  be  Justified. 

(2)  The  material  in  the  shell  is  also  considered  to  be  homogeneous, 
isotropic,  and  linearly  elastic.  Generally  speaking,  this  assumption  is  valid 
for  values  of  stress  below  the  so-called  proportional  limit  of  materials 
commonly  used.  In  addition,  the  thickness  of  the  shell  relative  to  its  radius 
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is  assumed  small;  this  permits  expression  of  all  stress  components  of  the 
shell  in  terms  of  a  function  which  describes  the  deflection  of  its  middle 
surface.  This  deflection  must  satisfy  a  linear  partial  differential  equation 
with  the  appropriate  boundary  conditions. 

(3)  The  incident  stress  wave  considered  is  either  a  plane  dilata- 
tional  or  a  plane  distortional  (shear)  wave.  Under  actual  conditions,  both 
waves  are  propagated  with  the  shear  lagging  the  dilatational  wave.  The 
combined  effect  for  elastic  conditions  may  be  determined  through  the  principle 
of  superposition. 

(4)  The  radial  and  tangential  particle  velocities  of  the  medium 
at  the  boundary  are  equal  to  that  of  the  shell.  This  is  the  continuity 
relation  insuring  that  the  shell  and  the  medium  are  in  contact  with  no 
relative  slip  occurring  at  the  boundary. 

(5)  Any  additional  mass  within  the  shell  is  assumed  to  be  dis¬ 
tributed  symmetrically  about  the  axis.  The  significance  of  this  assumption 
is  found  in  the  development  of  the  equations  of  motion  to  account  for  any 
additional  mass  located  within  the  shell. 

Other  assumptions  a  re  presented  in  the  formal  development  of  the 
mathematical  expressions  used  to  describe  the  behavior  of  shell  and  medium. 

1.4  Method  of  Approach 

The  two  components  of  shell  displacement,  radial  and  tangential, 
are  written  in  terms  of  Fourier  series  from  which  expressions  giving  the 
strain  energy  and  kinetic  energy  of  the  shell  in  vacuo  are  derived.  The 
equations  of  motion  are  deiived  from  Lagrange’s  equations  in  terms  of  the 
displacement  functions  and  forces  acting  on  the  shell.  Forces  on  the  shell 
result  from  the  stresses  in  the  medium  at  the  boundary.  Stresses  in  the 
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medium  are  taken  to  be  the  sum  of  the  stresses  due  to  the  incoming  stress  wave 
expressed  in  terms  of  Fourier  series  whose  coefficients  are  known,  and  those 
due  to  the  reflected  and  diffracted  effects  expressed  in  terms  of  a  pair  of 
displacement  potentials  representing  waves  diverging  from  the  axis  of  the 
shell.  The  form  of  these  potentials  as  derived  by  Lamb  (4)  ,  is  in  terms  of 
sine  and  cosine  series,  the  nature  and  treatment  of  which  has  been  studied 
by  Paul  (6) - 

The  equations  of  motion,  described  as  a  pair  of  coupled  integro- 
differential  equations,  are  solved  modewise  using  a  numerical  technique 
known  as  the  Newmark  Beta  Method  (5)  which  permits  determination  of  the 
coefficients  of  the  potentials  and  values  of  acceleration,  velocity,  and 
displacement  of  the  shell. 

The  solution  obtained  is  compared  with  values  obtained  from  a  series 
expansion  of  the  equations,  which  is  valid  for  short  times,  and  with  the  static 
solution  based  on  the  theory  of  elasticity  which  the  machine  solution  should 
approach  asymptotically. 

1.5  Previous  Work 

The  problem  stated  above  has  been  the  subject  of  a  recent  report 
by  Baron  (l),  whose  analysis  consists  of  first  solving  for  the  displacements 
caused  at  the  boundary  of  an  unlined  cylindrical  cavity  subjected  to  a  plane 
stress  wave.  This  is  done  through  an  integral  transform  approach,  the  solu¬ 
tion  of  the  transformed  equations  being  expressed  using  Hankel  functions. 

The  evaluation  of  the  inverse  transform  is  accomplished  only  with  great 
computational  effort.  Values  of  displacements  obtained  are  then  used  as 
influence  coefficients  in  determining  the  displacements  of  the  shell. 


*  Numbers  in  parenthesis  refer  to  the  corresponding  entry  in  the  Bibliography. 
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Solutions  obtained  for  two  different  sets  of  parameters  by  Baron  are  compared 
in  Section  4.  6  with  results  obtained  by  the  method  of  solution  outlined  in 
this  report. 

The  study  by  Paul  (6)  consisted  of  analyzing  the  effect  of  a  plane 
stress  wave  incident  on  an  unlined  cylindrical  cavity  in  an  elastic  medium. 
The  reflected  and  diffracted  waves  are  described  in  terms  of  displacement 
potentials  which  represent  outgoing  shear  and  dilatational  waves.  A  method 
was  developed  for  determining  values  of  these  potentials,  and  a  similar 
method  is  used  in  this  report. 


1.6  Notation 


Notation  is  defined  throughout  the  text  where  it  first  appears; 
however,  the  following  list  summarizes  the  main  uses  of  certain  symbols.  In 
discussions  of  special  topics  other  meanings  may  be  ascribed  to  the  symbols, 
at  which  time  they  will  be  redefined. 


A  =  Area  of  cross  section  of  the  shell,  per  unit  length. 

A  .  B  .  C  =  Coefficients  of  Fourier  series  for  stresses  in  the 
n'  n  n  ,, 

medium. 


A 


n 


y 


C 


n 


Coefficients  of  Fourier  series  for  total  stresses 
in  the  medium. 


"V 


at 


m 


Weighting  factors. 


a 


n' 


b 


n 


Generalized  coordinates  for  the  displacements  of 
the  medium  at  the  boundary. 


a 


ns' 


b 


ns 


Generalized  coordinates  for  displacements  of  the 
shell. 


c 


1' 


Velocities  of  wave  propagation,  dilatational  and 
shear,  respectively . 


d 


Distance  from  neutral  axis  of  shell  to  its  outermost 
fiber. 
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E,  E  =  Moduli  of  elasticity  for  the  medium  and  shell, 

*  respectively.  Bar  over  the  symbol  refers  to  the 

plane  strain  modulus. 

e  =  Volumetric  strain. 


F,  G 


I 

k 


m 

m' 


=  Generalized  coordinates  of  the  displacement 
potentials . 

=  Moment  of  inertia  of  the  shell,  per  unit  length. 

=  Parameter  which  relates  to  the  shape  of  a  time 
dependent  stress  wave. 

=  R/r,  ratio  of  radii. 

=  ratio  of  velocities. 

=  Mass  of  the  shell  per  unit  surface  area. 

»  Additional  mass  within  the  shell. 


n 


Mode  number. 


Generalized  forces. 


R 


Radius  of  the  shell. 


V  Sn 

r,  0 

T 

Tf 

t 

U 

u 

U,  V 


u 


=  Multiplying  factors. 

=  Polar  coordinates . 

=  Kinetic  energy  of  the  shell. 

=  Kinetic  energy  of  additional  mass  within  the  shell. 
=  Thickness  of  the  shell;  also,  time. 

=  Strain  energy  of  the  shell. 

=  Displacement  vector. 

=  Components  of  the  displacement  vector  in  the  radial 
and  tangential  directions,  respectively. 

=  Radial  and  tangential  displacement  components  of 
the  shell. 


V  uy 


Components  of  the  displacement  vector  in  the  x  and  y 
directions ,  respectively . 


x,  y,  z 


Rectangular  coordinates. 
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a  E 
sn  s 

^  ,  non-dimens  ionalized  radial  component  of 

p  displacement. 

b  E 
sn  s 

^  — —  ,  non-dimensionalized  tangential  component 
p  of  displacement. 

Strains  along  the  x  and  y  axes,  respectively. 
Circumferential  strain. 

A  variable  of  integration. 


,  a  parameter, 
a  parameter. 

,  a  parameter. 

a  parameter. 

,  a  parameter. 


,  a  parameter. 

Polar  coordinates. 

Position  angle  of  the  incident  wave. 

Curvature  of  the  shell. 

Lame4  constants. 

Poisson's  ratio. 
v 

1-v 

Arguments  of  the  functions  F  and  G,  respectively. 

Mass  density  of  the  medium. 

Mass  density  of  the  shell. 

Normal  components  of  stress  parallel  to  the 
x  and  y  axis,  respectively. 


Shearing  stress  component. 
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rr’  °99 

Radial  and  tangential  normal  stresses  in  polar 
coordinates.  Additional  subscripts  s  and  m, 
when  used,  refer  to  the  shell  and  medium. 

't9 

= 

Shearing  stress  in  polar  coordinates. 

'sb 

= 

Bending  stress  in  the  shell. 

r 

smax 

- 

Maximum  stress  in  the  shell. 

r 

P 

= 

Amplitude  of  incident  dilatational  wave. 

r 

s 

= 

Amplitude  of  incident  shear  wave. 

T 

= 

Non-dimensionalized  time  in  the  case  of  the 
incident  dilatational  wave. 

= 

Non-dimensionalized  time  in  the  case  of  the 
incident  shear  wave. 

>,  * 

= 

Displacement  potentials,  dilatational  and  shear 
respectively. 
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CHAPTER  IT 
BASIC  EQUATIONS 


2,1  Equations  for  the  Medium 

Tht  differential  equation  of  motion  of  a  particle  in  a  homogeneous, 
isotropic,  and  linearly -elastic  medium  in  terms  of  its  displacement  vector  u 
is  given  by  Kolsky  (j)  in  the  form 


(\  +  2u)  V  V  *  u  -  \i V  x  (7  x  u) 


(2-1) 


vhere  A  and  u  are  the  Lame'  constants,  and  p  is  the  density  of  the  medium.  The 
vector  u  may  be  expressed  as  the  sum  of  two  displacements,  the  gradient  of 
a  scalar  potential  and  the  curl  of  a  vector  potential 


u  »  79  +  If 


(2-2) 


Here  9  is  a  potential  giving  rise  to  an  irrotational  displacement  and  ¥  a 
potential  leading  to  an  equivoluminal  displacement. 

If  the  wave  equations  (2-3)  are  satisfied  and  u  expressed  as  in 
(2-2),  the  equations  of  motion  (2-1)  are  automatically  satisfied. 


2  « 
57  9 


1  a  q> 
c2  dt2 


4  * 


2  i*2 

ot 


(2-3) 


where  the  velocities  of  wave  propagation  are 


(2-4) 


for  the  dllatational  and  shear  waves,  respectively. 
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Since  the  problem  is  essentially  a  two-dimensional  one  wherein  only  the 
rotation  about  the  cylinder  axis  is  considered,  the  expression  for  the 
propagation  of  shear  waves  can  be  written  in  terms  of  a  scalar  potential 
function 


2 

V  * 


_1_  cQ 

"  2  x+2 

c2  dt 


with  \|r  =  V 


(2-3a) 


cp  and  are  functions  of  x,  y,  and  t.  The  components  of  the  displacement 
vector  u  can  be  expressed  as 


u  & 

x  ox  ay 


(2-5) 


U  = 

y 


Stress  components  are 


a  =  \e  +  2pe 

XX  X 

a  =  UT 
xy  xy 

o  =  \e  +  2pe 

yy  y 


(2-6) 


where  e  represents  volumetric  strain,  €  and  e  strain  along  the  x  and  y  axis 

x  y 

respectively,  and  x  the  shearing  strain.  The  strain  components  for  small 
xy 


displacement  are 


du 

_ X 

sr 


du 

‘  ^ 


(2-6a) 


xy 


du  du 

5s ’S1 


In  terms  of  the  potential  functions,  the  equations  for  the  stress  components 


become 
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o  =  (i 
xy 


(2  it. .  ii  +  £i) 

to*  dx2 


In  polar  coordinates 

2  a2  i  a  i  a2 

7  'dr2  ?5?  r2  d#2 


and  the  wave  equations  expressed  in  polar  coordinates  are 


£2.  +  i&.  +  ±£i  =  ±is> 

>  2  r  dr  2  .fl2  2-^2 

or  r  00  dt 

a2t  i  at  i  a2*  1  a2t 

•v  2  r  Sr  2  -^2  2  v.2 

or  r  0©  c2  ot 


Displacements  in  polar  coordinates  become 

1  d4 


U=£  + 


|§ 

r  do  or 


(2-7) 


(2-8) 


(2-9) 


where  u  and  v  represent  radial  and  tangential  components  of  displacement, 
respectively. 


Stress  components  may  be  written  as 

SQ  "  , 

r 


v  -  (d  +  2n>[lf  +  risk '  7  H]  *  7  [if  +  r  F  H  •  H§] 


r0 


90 


(2-10) 
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II 

II 

II 

II 

II 

ill 

III 

II 


( 
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where  a^}  o^q,  a^,  represent  respectively  the  radial,  shear  and  circum¬ 
ferential  (hoop)  stresses  in  the  medium. 

In  the  following  development  of  the  equations,  the  incoming  wave, 
considered  a  step  function  of  time,  is  either  a  dilatational  or  shear  wave. 

The  step  wave  was  chosen  as  a  convenient  form  from  which,  through  the  principle 
of  superposition  using  Duhamel’s  integral,  the  effects  of  any  time  dependent 
wave  may  be  approximated.  It  is  shown  by  Kolsky  (3)  that  for  either  the 
incoming  dilatational  or  shear  wave,  the  reflected  wave  from  a  plane  body  must 
have  both  dilatational  and  shear  components  in  order  to  satisfy  boundary  con¬ 
ditions.  The  same  holds  in  the  present  case. 

2.11  Dilatational  Wave 

The  orientation  of  the  incoming  dilatational  (compressional)  wave 
with  respect  to  the  spatial  coordinates  and  time  are  shown  in  Fig.  1.  Owing 
to  symmetry  about  the  x  axis,  the  radial  stress,  hoop  stress,  and  the  radial 
component  of  displacement  are  even  functions  of  0,  while  the  shear  stress  and 
tangential  component  of  displacement  are  odd  functions  of  0. 

In  the  expressions  previously  given,  the  dilatational  potential  may 
now  be  separated  into  two  parts 

*  =  V  +  *out 

where 

v  *  -  °P  ^  (x+cit} 

represents  the  potential  of  the  incoming  step  wave,  and  <p  represents  the 
potential  of  the  outgoing  dilatational  wave. 

The  displacement  and  velocity  of  a  particle  behind  the  plane  wave 
front  can  be  derived  from  the  potential  of  the  incoming  wave,  with  notation 
defined  in  Fig.  1 
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ux  *  “  ap  ^ ' (x+cit^ 


a 


The  velocities  in  polar  coordinates  become 

o 


u  =  -  JP  cos  9 
pci 

0 

v  =  —2—  sin  9 

QZ1 


(2-11) 


(2-12) 


The  radial  and  tangential  displacements  at  the  boundary  r  =  R  due  to  the 
incoming  plane  stress  wave  are  expanded  in  Fourier  series  as  a  function  of 
time  as 


u(R,t,e)  =  £  an(R,t)  cos  n0 
n=0 

(2-15) 

00 

v(R,t,6)  =  £  bn(R,t)  sin  n9 
n=l 

and  the  velocities  as 


u(R,t,e)  =  Y  cos  nS 

n=0 

oo  (2-14) 

v(R,t,9)  =  £  bn(R,t)  sin  n0 
n=l 

The  Fourier  coefficients  are  determined  as  a  function  of  a  measure  of  the 
degree  of  envelopment  of  the  shell  by  the  incoming  stress  wave  (Fig.  l),  by 
evaluation  of  the  integrals 
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cos  0  cos  n0  d0 


sin  0  sin  n 9  d0 


The  indicated  integrations  are  performed  using  the  orthogonality  properties 
of  the  sine  and  cosine  functions  to  get  the  following  result 


*n(in) 


a 

—L. 

«pc. 


3n(in) 


a 

_J2_ 

*PCX 


sin  9 


1 

sin  20, 


*1  + 


sin(n-l)0^  sin(n+l)01 


n-1 


n+1 


sin  20, 


sin(n-l)01  sin(n+l)01 


n-1 


n+1 


n  =  0 


n  =  1 


n  —  2 ,  3>  • • • 
(2-15) 

n  =  1 

n  *  2,  3>  •  •  • 


The  subscript  (in)  has  been  added  to  designate  the  coefficients  due  to  the 
incident  wave*  0^  is  a  time  dependent  variable  which  varies  from  0  to  it 
during  the  transit  of  the  incoming  wave  across  the  cavity,  after  which  it 
becomes  a  constant  equal  to  *• 

The  stresses  in  the  medium  behind  the  wave  front  are  as  shown  in 
_  * 

Fig.  1,  where  v>  a  ratio  relating  the  stress  in  a  direction  perpendicular  to 
the  direction  of  wave  propagation,  is  derived  from  the  assumption  that  there 
is  no  strain  behind  the  wave  front  in  the  direction  parallel  to  the  wave  front. 
It  is  related  to  Poisson* s  ratio,  v,  of  the  medium  in  the  following  way 


*  The  bar  over  the  symbol  does  not  indicate  a  vector  quantity. 
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v  = 


1  -  V 

The  stresses  behind  the  wave  front  are  in  polar  coordinates 

ffrr  =  "  ap(cos20  +  v  sin2e) 

a  a  =  a  (^5^)  sin  20 
rd  p  x  2  ' 

2^2 

oQq  =  -a^(sin  6  +  v  cos  0) 


(2-16) 


(2-17) 


Stresses  in  the  medium  as  a  function  of  radius  and  time  are  now 
expressed  in  terms  of  Fourier  s tries 


’rr^'®)  =  I  An(r’t}  cos  n$ 


n=0 


so 

(r,t,6)  =  Y.  Bn(r't)  sin  n0 


(2-18) 


n-1 


a66  (r,t,0)  =  £  Cn(r,t)  cos  n 9 


n=0 

and  the  coefficients  determined  in  the  same  manner  as  used  previously  for 
velocities  are 

f  .  .1.77. 

n  =  0 


(l+v)ei  +  (^)  sin  2B1 


An(in)  =  "  2«  *  1  (1-v)®i  +  (1+v)  sin  20i  +  (^)  8in  n  =  2 


(l-v)sin(n-2)01  (l-v)sin(n+2)01  2(l+v)  sin  n©1 

— -  +  - 


(2-19) 


n+2 


n 


B 


.  i.(i-7) 


sin  40 

ei  r-2 


h(in)  ~  2* 


sin(n-2)0^  sin(n+2)0^ 


n-2 


n+2 


n  =  1,3, ^,5,  • 


n  =  2 


n  =  1,3, 4, 5,. 
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(l+v)^  -  (ig*)  sin  29 1  n  =  0 

•  1  (v-l)61  +  (l+v)  sin  20^  +  (^jp)  sin  4^  n  =  2 

(v-l)sin(n-2)01  (v-l)sin(n+2)0^  2(l+v)sin  n©1 

n^2  +  n+2  +  n 

n  =  1, 3;  5 >  •  •  • 

\ 

where  the  subscript  (in)  again  refers  to  the  incident  stress  wave. 

As  noted  earlier,  the  reflected  and  diffracted  wave  is  represented 
as  a  diverging  or  outgoing  cylindrical  wave  from  a  line  source  whose  origin  is 
the  axis  cf  the  cylinder,  and  includes  both  dilatational  and  shear  components. 
The  outgoing  potential  functions  are  expressed  as  infinite  series,  each  term 
of  which  will  henceforth  be  called  a  mode 


00 

*out  =  Z  fn(r’t}  C0S  n0 


n=0 


out 


00 

"  Z 


sin  n© 


(2-20) 


n=l 


It  is  shown  in  the  appendix  that  the  coefficients  are  of  the  form 

00 

f  =  /  ^(t  -  cosh  u. )  cosh  nu,  du, 

n  cn  JQ  ex  1'  11 

00 

«n  =  £  °n(t  -  t;  cosh  u2)  C08h  nu2  du2 


(2-21) 


C2 


where  F°(t  -  cosh  u  )  and  Gn(t  -  —  cosh  u  )  represent  the  nth  order 

C1  X  c2  2 

derivative  of  the  respective  functions,  n  being  equal  to  the  mode  number 


considered. 
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To  express  the  velocity  components  of  the  medium  at  the  boundary 
due  to  the  outgoing  waves  as  in  equations  (2-14),  one  differentiates 
equations  (2-9)  with  respect  to  time  and  substitutes  the  expressions  given 
above  for  the  potential  functions.  This  results  in  the  following  for  the 
velocity  coefficients 


00 

4n(°ut)  ’  "  L  f 2<‘l>  cosh  U1  cosh  nul  dul 

C1 

00 

+  J  Gn+2(l2)  sinh  u2  sinh  nu 2  dUg 

c2  0 


n(out) 


00 

=  -  J  Fn+2(61)  sinh  sinh  nu1  dux 


00 

^1"  1  Gn+2(62}  C0Sh  u2  COSh  nu2  dU2 


c“  w0 


(2-22) 


where  (61)  and  (Sg)  represent  the  arguments  of  the  respective  functions. 
The  subscript  (out)  refers  to  the  outgoing  wave. 

The  stresses  in  the  medium  due  to  the  outgoing  waves  are  also 
expressed  in  series  form  as  in  equations  (2-18),  the  coefficients  being 
found  by  substitution  into  equations  (2-10) 


<-Unn 

^n(out)  ”  n+2 
C1 


l  ^2“i> 


\  2 

cosh  nu..[—  +  2  cosh  u,  ]  du.. 

l‘(i  11 


0ntS({2) 


sinh  2u2  sinh  nu2  dUg 


(2-23) 


18 


Bn(out) 


sinh  2u^  sinh  nu^  du^ 


L±Zjl 

n+2 


00 

l 


0n+2(i2) 


cosh  2u^  cosh  nu2  du^ 


C 


n(out) 


00 

f  Fn+2(Si)  cosh  nui  f“  "  2  sinh2U;L]  d^ 
C1  0 

00 

+  J  Gn+2U2)  sinh  2u2  sinh  nu£  du2 

C2  ° 


2.12  Shear  Wave 

The  orientation  of  the  incoming  plane  shear  wave  with  respect  to 
the  spatial  coordinates  and  time  are  shown  in  Fig.  2.  In  this  case  the  shear 
stress  and  tangential  component  of  displacement  must  be  even  functions  of  9 
while  the  radial  stress,  hoop  stress,  and  radial  component  of  displacement 
are  odd  functions  of  9 . 

The  incoming  wave  expressed  in  terms  of  a  shear  potential  function 


is 


♦  .  =  o  ")C  (x  +  c0t) 

in  s  s'  2 

from  which  the  velocity  of  a  particle  behind  the  wave  front  in  polar 
coordinates  may  be  derived  as 
a 


(2-24) 


u  = 


s 

pc. 


V  = 


pc. 


sin  q 


cos  9 


(2-25) 


The  stresses  in  the  medium  behind  the  wave  front  are  as  shown  in  Fig.  2.  In 
terms  of  polar  coordinates  they  are 
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a  =  a  sin  20 
rr  s 

ar9  “  as  008  20  (2-26) 

aee  “  “a8  8in  20 

Velocities  and  stresses  due  to  the  incoming  plane  stress  wave  as  functions 
of  radial  distance  and  time  are  expanded  in  terms  of  Fourier  series  in  a 
manner  similar  to  that  of  the  dilatational  wave,  except  that  because  of  the 
difference  in  symmetry,  the  sine  and  cosine  terms  are  interchanged.  The 
coefficients  are  also  found  in  a  similar  manner  and  are  given  here.  The 
coefficients  for  the  velocity  series  are 


Si(in)  ”  xpc-  x  \ 


sin  20, 


*1" 


sinCn-l^  sin(n+l)0^ 


n-1 


n+1 


n  =  1 


n  =  2,3,... 

(2-27) 


3n(in)  itpcg  x 


sin  0^ 

sin  2 &1 


sin(n-l)0^  sin(n+l)0^ 
n-1  +  n+1 


n  =  0 

n  =  1 


n  =  2,3,... 


and  the  coefficients  for  the  stress  series  are 


a 

A  /.  v  =  -r  x  < 
n(in)  «  N 


sin  40 

sin(n-2)0^  sin(n+2)0^ 


n-2 


n+2 


B 


sin  20, 


sin  40, 


n(in) 


-2.  X  \  Q  + 
n  1  ’ 


sin(n-2)01  sin(n+2)0;L 


n-2 


n+2 


n  =  2 

n  =  1,3,4, 3,.*. 


n  »  0 

n  =  2 

(2-28) 

n  —  1,3, 4, 5, . . . 
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Cn(in)  ”  '  ^n(in) 


0^  is  a  time  dependent  variable  similar  to  that  for  the  dilatational  wave 
except  that  now  it  is  a  function  of  the  position  of  the  incoming  shear  wave, 
which  travels  with  velocity  Cg. 

The  outgoing  potential  functions  are  expressed  modewise  as 


%t  =  I  8in 


nQ 


n=l 


out 


w 

“  X  8n(r>t) 


cos  n 0 


n=0 


(2-29) 


where  the  coefficients  are  of  the  form 


n  °° 

f  =  jp-  ■  /  F^(  t  -  —  cosh  u, )  cosh  nu,  du. 

n  cn  JQ  cx  1  11 

n  00 

gn  =  J  Gn(t  -  -2-  cosh  u2)  cosh  nug  dug 


(2-30) 


The  terms  have  the  same  significance  as  in  the  case  of  the  dilatational  wave. 
Modal  coefficients  are  derived  as  "before.  The  coefficients  for  the  velocity 
terms  are 

00 

Vout) =  -  L  ,p+2<fti)  cosh  ui cosh  nUl  dUl 


00 

^+1  X  °n+2(52)  8inh  U2  8inh  nJ2  du2 


c2  -  -0 


(2-31) 
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00 

Vout)  ■  j[  sinh  ui sinh  nui  dui 

C1  u 


00 

^Jl  X°n+2(62)  C08h  u2  C08h  nu2  du2 
c2  ° 


and  the  coefficients  for  the  stress  terms  are 


00 

An(o»t)  '  ^  L  F”"2(81)  c08h  nul^  *  2  C08h2  "l1  aul 


C1  0 


CO 

^'~n+2^  J  °n+2^2^  Sinh  ^2  Slnh  nu2  du2 


V  -0 


B 


n(out) 


00 

=  ■  X  Fn+2(si)  sinh  sinh  nui  dui 


c1_  w  0 


n  00 

’  ^~n+2^  J  °n+2^2^  COSh  ^2  COsh  nu2  du£ 


c2  "  0 


(2-32) 


00 

Cn(out)  ’  ^  L  c08h  nultj  -  2  8lnh2  V  a“l 

C1  ^ 

00 

-  ^ ~n+2^  J  °n+2^2^  8lnh  2u2  8inh  nu2  du2 


c2  '  0 


2.2  Equations  for  the  Shell 

In  the  following  derivations  of  the  equations  of  motion  for  the 
shell,  it  is  assumed  that  the  thickness  of  the  shell  is  small  in  relation  to 
its  radius.  This  assumption  permits  description  of  the  behavior  of  the  shell 
in  terms  of  its  middle  surface,  and  effectively  concentrates  all  of  the  mass 
in  a  line  thickness.  This  should  be  kept  in  mind  when  interpreting  the  results. 
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The  effects  of  shear  and  rotatory  inertia  are  not  considered,  since  they 

affect  only  the  very  high  modes,  which  are  not  considered  in  this  study. 

The  deformation  of  any  point  (0)  on  the  shell  can  then  be  described 

by  specifying  the  components  of  displacement  (Pig.  3) - 

u  ■  radial  displacement,  positive  outwards 
s 

v  =  tangential  displacement,  positive  in  the  direction 
of  increasing  0. 

The  strain  energy  of  the  shell  per  unit  length  can  be  expressed  as 
=  T  2*  =  .  2x 

u  aJrf0  (lO^de  +  -§r  (e0)2 Rde  (2-33) 


where 


Es  =  "plane  strain"  modulus  of  elasticity  of  the  shell 
I  =  moment  of  inertia  per  unit  length  of  shell 
A  =  area  of  the  cross  section  per  unit  length  of  shell 
k  =  change  in  curvature  of  the  cross  section 
=  extensional  strain  in  the  0  direction 
R  =  radius  of  the  cylinder 


Flugge  (2)  gives  the  expressions  for  the  curvature  and  extensional  strain 
terms  as 

,  d2u 

k  =  ^2  +  us) 

R2  b92  s 

and 

H  ‘  i(5T*  us> 

correct  to  first  order  terms.  The  strain  energy  equation  in  terms  of  the 
components  of  displacement  then  becomes 

2n 


2  n  .2 
o  u 


Ei  p  o 

U  =  A  /  (- 

2R-7  J  0  b9 


\2 

2  +  Us}  +  2R 


P  bv  2 

(sr +  us)  de 


(2-3^) 
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The  kinetic  energy  of  the  shell  can  be  expressed  as 
2ff 

T  •  1  /0  *  *s>  M9 

where 

m  =  mass  of  shell  per  unit  surface  area 


(2-55) 


2.21  Dilatational  Wave 

In  view  of  the  symmetries  of  the  problem,  the  displacement  com¬ 
ponents  of  the  shell  can  be  expanded  in  terms  of  the  following  Fourier  series 


us(e,t)  =  £  asn(t)  cos  nO 
n=0 

(2-56) 

00 

v8(e,t)  =  £  bsn(t)  sin  n9 
n=l 


Equations  for  the  strain  and  kinetic  energies  of  the  shell  are  then  expressed 
as  the  following  quadratic  forms  in  the  generalized  displacements  and 
velocities 


E  Iff  E  Ax 


u  =  (-V  + 

R-' 


E  Aff  _  E  I#  P  „  0  0  EA*  n 

~ — )  a^  +  — V  a^  (l-n^)  +  y  (nb  +  a  ) 

R  '  so  ^5  Z_i  an'  2R  sn  sn' 


o  E  Aff 

2  +  — 
2R 


n=l 


n=l 


(2-37) 


_  ^ .  2  nutR 

T  «  mnRa  + 

so  2 


jv^> 

sn  sn' 


n=l 


The  Lagrangian  equations  of  motion  in  terms  of  the  generalized  displacement 

coefficients  a  and  b  are 
sn  sn 


d  /dT 
dt  '5a 

sn 


) 


Q 

ST-  =  Si 

sn 


d  /ftr 

dt  (5T~ 
sn 


) 


dU 

5T" 

sn 


Si 


(2-58) 
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where  and  ^  are  generalized  forces  corresponding  to  the  displacement 
terms  agn  and  bgn,  respectively.  The  following  equations  of  motion  result 


El  E  A  Q 

••  /  8  S\  0 


El  _  2  BA  Q 

a  +  "  i  (l-n  )  a  +  —?  (a  +  rib  )  =  _ 

sn  ^4  '  '  sn  ^2  v  sn  sn'  m«R 


nE  A  o 

b  +  — i—  (a  +  nb  )  =  — 

sn  ^2  '  sn  sn'  m«R 


n  =  0 

(2-39) 

n  =  1,  2,  ... 


n  =  1,  2, 


Now  to  determine  the  generalized  forces  and  the  principle 

of  virtual  work  is  applied.  Consider  a  virtual  displacement  corresponding 

to  an  increment  8a  of  coordinate  a 

sn  sn 

Virtual  work,  by  definition  =  •  8agn 

The  external  forces  acting  on  the  shell  are 


00 

•IX 


or  =  )  A  cos  n0 

rr 

n=0 
00 


r« 


sin  n0 


n=l 


where  and  B  represent  the  sum  of  the  forces  on  the  shell  due  to  the 
n  n 

incoming  and  outgoing  stress  waves 


An  =  An(in)  +  An(out) 
Bn  =  Bn(in)  +  Bn(out) 


The  coefficients  An  and  Bn  have  been  discussed  previously, 
a  virtual  change  of  the  generalized  coordinate  agn  is 

2n 


Virtual  Work  = 


—  2 

A  cos  n0  •  6a 
n  sn 


•  Rd0 


The  work  due  to 
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The  generalized  force  term  is  determined  hy  equating  the  two  expressions  of 
work,  from  which  we  get 

Q  =  2«R  A  n  =  0 

o  o 

=  >tR  An  n  =  1,  2,  ... 

and  by  a  similar  computation 

*n  ”  *R  \ 

The  equations  of  motion  now  become 


El 


E  A 


8  +  *i  (l-n  )  a  +  (a  +  nb  )  =  - ~ 

sn  ^4  '  '  sn  ^2  '  sn  sn'  m 


nE  A  B 

B  (a  +  nb  )  = 

sn  _2  sn  sn'  *■ 
tnn 


m 


(2-40) 


2.22  Shear  Wave 

Because  of  the  difference  in  symmetry  associated  with  the  shear 
wave  as  compared  to  the  dilatational  wave,  displacement  components  in  the 
case  of  the  Incoming  shear  wave  are  written  as 


sin  nB 

n=l 

00 

v8(»,t)  =  £  b8n(t)  cos  n 0 
n=0 


(2-41) 


Expressions  for  the  strain  and  kinetic  energies  in  terms  of  the  generalized 

coordinates  a  and  b  now  become 
sn  sn 


u  = 


a  ) 
sn' 


2 


(2-42) 
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_  -2  mnR 

T  =  mflHb  +  —t— 
so  2 


00 


•2 

a  +  b  ) 
sn  sn' 


n=sl 

The  equations  of  motion  determined  in  a  manner  similar  to  that  of  the 
dilatational  wave  are  then 
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2.23  Effect  of  Additional  Mass 

The  equations  of  motion  derived  thus  far  assumed  no  mass  within  the 
shell.  The  simplest  way  to  consider  the  effect  of  additional  mass  is  to 
assume  that  it  is  distributed  within  the  shell  symmetrically  about  the  axis. 
The  total  quantity  can  be  assumed  equal  to  2JtRm* ,  with  m1  the  magnitude  of 
the  added  mass  measured  in  terms  of  the  surface  area  of  the  shell.  If  it  is 
assumed  to  move  with  the  same  velocity  as  the  mass  center  of  the  shell,  its 
kinetic  energy  can  be  expressed  as 

2%  2 

T*  =  I  J'  (-\ig  cos  0  +  vg  sin  Q) d$j  (2-WO 

for  the  case  of  the  incoming  dilatational  wave.  If  u  and  v  are  replaced  by 

s  s 

their  Fourier  expansions,  and  if  the  indicated  integration  is  carried  out, 
the  kinetic  energy  T1  due  to  the  additional  mass  may  be  represented  as 


T'  = 


m '  JtR  / .  2 


(a 


si 


-  2a  ,b 
si  si 


This  kinetic  energy  uiust  be  added  to  the  kinetic  energy  of  the  shell  which 
was  derived  earlier.  Since  the  above  expression  affects  only  the  n  =  1  mode, 
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the  other  modes  need  not  be  considered.  With  this  additional  kinetic  energy, 
the  equations  of  motion  for  the  n  =  1  mode  become 


,.  m’  /..  r  \  V  /  \  A1 
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If  now  we  add  the  above  equations,  we  get 
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By  making  the  appropriate  substitutions  for  agl  and  b  the  equations  of 
motion  can  be  written  as 


E  A  « 
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sl  -2  v  sl  sl7  m 
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mR 


.  .  .  m'  1  t  2m-m’ 

,2  '  sl  sl'  m  ’  2(m-Hn' )  m  2(m+m' ) 


Note  that  these  equations  are  similar  in  form  to  the  equations  derived  earlier 
without  the  additional  mass;  the  effect  of  the  additional  mass  merely  alters 
the  right-hand  side. 

In  the  case  of  the  incoming  shear  wave,  the  kinetic  energy  of  the 


additional  mass  is 
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2ji  2 

T'  =  ^  fQ  ^8  8ln  9  +  ^8  COB  9)  d9]  (2-Vf) 

and  by  a  similar  method,  the  equations  of  motion  for  the  n  =  1  mode  are 
written 

..  V  /  .  \  =tl  .  2m+m'  _  \  _  m’ 
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The  significance  of  the  additional  mass  on  the  numerical  results  is  discussed 
in  a  later  chapter.  Inclusion  of  a  flexible  support  for  the  additional  mass 
is  also  possible  by  changing  the  above  equations;  however,  this  problem  is 
not  treated  here. 

2.3  Boundary  Conditions 

It  is  assumed  that  the  shell  is  attached  to  the  medium  permitting 
no  differential  displacements  at  the  boundary  between  the  two.  Thus,  con¬ 
tinuity  of  stresses  and  displacements  are  maintained  at  the  boundary.  A 
convenient  way  of  satisfying  the  boundary  condition  of  equal  displacements 
of  the  shell  and  medium  is  to  equate  the  corresponding  velocity  components 
of  the  shell  and  medium.  This  is  done  modewise  for  both  the  dilatational 
and  shear  incoming  waves  by  merely  equating  the  coefficients  of  the  velocity 
terms  as  follows 

*sn  "  *n(in)  +  *n(out) 

V  =  \(in)  +  *n(out> 


(2-49) 
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2.k  Summary  of  Equations  in  Non-Dlmensionalized  Form 

It  is  convenient  to  express  the  equations  derived  earlier  in 
non-dimens ionalized  form  by  introducing  certain  dimensionless  variables  and 
parameters .  Let 


(2-50) 


be  the  non-dimensionalized  form  of  the  displacement  and  time  variables. 
Parameters  are  expressed  as 


\ 


\  P 

«  _  (l+v)(l-2v) 
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(2-51) 
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The  last  two  parameters  given  above  become  unnecessary  when  considering  an 
unstiffened  shell  since  for  this  case  A  and  I  are  given  by 


and 


i  _  _L  t5 

1  12  1 


(2-52) 


Although  the  equations  were  derived  considering  an  unstiffened  shell  of 
uniform  thickness,  the  response  of  a  shell  with  thin,  closely  spaced  stiffeners 
can  also  be  approximated  by  specifying  A  and  I  separately  from  t.  The 
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equations  of  interest  in  non-dimensionalized  form  are 

dn  "  P1n(V  +  ^n(P»0) 

Pn  -  +  ^2n(r'G) 

(2-53) 

“n  ■  *Vai d'»»>  *  PW  +  «n(F>0) 

K  ■  +  PV(91>  +  <IVP’0> 


where  FI  .  F2  .  F3  <  and  P4  are  functions  whose  values  are  determined 
n  n  n  n 

directly  from  the  position  of  the  incident  wave.  Qln,  Q2r,  Q3n>  and  Q4n  are 
functions  of  the  outgoing  shear  and  dilatational  waves;  and  Nlfl  and  N2r  are 
functions  of  the  displacements.  These  various  functions  in  non-dimensionalized 
form  are  derived  from  equations  given  earlier.  Given  here  are  those  for  the 


case  of  the  incoming  dilatational  stress  wave. 
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2.5  Other  Equations  of  Interest 


Dilatational  Wave 

The  hoop  stress  in  the  shell  is  determined  from  the  relation 


e 


e 


where 


dv 

<3T  +  u-} 


By  substitution  of  the  Fourier  expansion  of  the  displacement  components  the 
hoop  stress  can  be  expressed  modewise  as  a  fraction  of  the  absolute  amplitude 
of  the  incoming  stress  wave,  as 


To2!  =  X  [nPn  +  °n]  cos  n® 
p  n=0 


The  bending  stress  is  determined  from  the  relation 


E  d  d2u 


which  may  also  be  expressed  modewise  as 


(2-55) 


i-8t  =  d  ^  (n2-l)  an  cos  nfl 
P  n=0 

where  d  is  the  distance  from  the  neutral  axis 
fiber . 


(2-56) 

of  the  shell  to  its  extreme 


Stresses  within  the  medium  at  any  radius  r  may  be  determined  by 
summation  of  modal  stresses  as  follows 


cos  n0 


(2-57) 
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where  An,  Bn>  and  are  coefficients  previously  described  as  generalized 
stresses  of  the  incoming  and  outgoing  stress  waves. 


Shear  Wave 

The  corresponding  equations  for  the  incoming  shear  wave  are,  for 

the  shell 
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and  for  the  stresses 
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where  all  terms  are  as  defined  before. 


CHAPTER  III 


34 


METHOD  OP  SOLUTION 


5-1  General 

In  this  chapter  are  presented  the  numerical  techniques  used  in 
solving  the  problem.  Since  the  methods  employed  are  similar  in  theory  for 
both  the  incoming  dilatational  and  shear  waves,  only  the  solution  to  the 
incoming  dilatational  (compressional)  wave  will  be  discussed  in  detail  here. 
The  equations  (2-53)  to  be  solved  consist  of  two  pairs  of  coupled  integro- 
differential  equations  in  the  generalized  coordinates  of  the  shell  and  the 
potential  functions .  The  integral  terms  in  these  equations  contain  elements 
of  the  outgoing  wave  potentials,  both  dilatational  and  shear,  whose  values 
must  be  determined  through  application  of  the  boundary  equations  at  each 
instant  of  time  considered. 

In  the  numerical  solution,  time  is  taken  to  be  the  independent 
variable,  and  is  expressed  in  terms  of  the  half  transit  time  of  the  incoming 
wave  across  the  cavity,  t  =  tc^/R.  Increments  of  time  are  expressed  as  a 
fraction  of  the  half  transit  time;  if  N  is  defined  as  the  number  of  time 
steps  for  the  wave  to  travel  one  radius,  then  l/N  denotes  the  elaspsed  time 
after  the  ith  step  has  been  taken  (Fig.  4). 

The  position  angle,  can  be  defined  for  any  step  in  time  as 

=  arc  cos  (1  -  i/N)  (3-1) 

After  total  envelopment,  i.e.,  for  i/N  >  2,  0^  becomes  equal  to  ir. 

3*2  Numerical  Integration  of  the  Potential  Functions 

Lamb  (4)  notes  that  if  a  point  source  f(t)dz  is  located  at  the 
point  (x=0,  y=0,  z=z),  its  effect  at  a  distance  r  from  the  origin  in 
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the  xy  plane  can  he  represented  by  the  equation 
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<P 


f  2  2 

f (t .  'Ll ±2-)  dz 


4tt  *</  r®+z^  '"l 

which  may  be  represented  graphically  as  in  Fig.  5*  To  obtain  the  effect  of 
a  line  source  of  density  F(t)  on  the  z  axis,  spherically  symmetric  point 
sources  with  the  same  variation  in  time  can  be  assumed  to  be  situated  all 
along  the  axis,  and  their  combined  effect  represented  as  an  integral  over  all 
the  point  sources 


•J  r^+z^ 


”  ‘  £  I. T(t ' 


dz 


•J  r2+z^ 


(3-2) 


The  limits  of  minus  infinity  to  plus  infinity  are  shown  here  in  order  to 

represent  the  general  case.  It  will  be  shown  that  for  a  disturbance  beginning 

at  some  definite  time,  the  limits  of  the  integral  may  be  taken  as  finite. 

Now  let 

2  2  2/,  2,  2\ 

r  +  z  =r(l+z/r) 

2  2 

*  r  cosh  u, 

and  the  above  integral  can  be  written  as 

ap 

<P  -  /  p(t  -  cosh  u1)  du.  (3-3) 

J0  ci  11 

from  which  the  graphical  interpretation  in  Fig.  6  is  made. 

Consider  for  the  present  the  determination  of  the  integral  values 
at  the  boundary  r  =  R.  We  define  t  =  0  to  be  that  time  when  the  first  effects 
of  the  outgoing  stress  waves  (i.e.,  from  the  point  source  at  z  =  0)  reach  the 
boundary.  As  time  increases  the  effects  of  other  point  sources  progressively 
further  away  along  the  z  axis  will  reach  the  boundary  in  the  xy  plane.  Thus 
it  can  be  seen  that  though  the  limits  of  infinity  are  given  in  the  integrals, 
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for  finite  times  the  limits  can  be  represented  by  finite  values.  The  maximum 
value  of  the  variable  R  cosh  u^  in  equation  (3-3)  that  need  be  considered  can 
therefore  be  represented  as  a  function  of  time.  At  t  *  0 
R  cosh  u1m  =  R 

that  is,  u1m  is  zero  for  time  zero.  For  any  time  later 
R  cosh  u.|m  =  R  +  tc^ 


in  the  case  of  the  outgoing  dilatational  wave  (Fig.  6).  R  cosh  u1n|  locates 
the  point  farthest  along  the  z  axis  whose  wave  front  has  Just  reached  the 
boundary  in  the  xy  plane,  c^  being  the  velocity  of  wave  propagation.  In 
non-dimensionalized  form  this  can  be  written 

COSh  U]Lm  =  1  +  T  =  1  +  i 

where  the  terms  are  as  previously  defined.  Since  the  integration  is  over  all 
point  sources  whose  wave  front  has  reached  or  passed  the  boundary  in  the  xy 
plane,  the  integral  expression  of  equation  (3-3)  can  be  written  with  new 
limits 

ulm 

=  J  F(t  -  2-  cosh  u^)  du^ 

Let 

cosh  u,  *  1  +  {  (3-^) 


where  £  now  becomes  the  variable  of  integration  varying  in  value  from 
0  to  t.  In  terms  of  C  the  integnl  expression  becomes 

t 


<P  - 


r 

> I l+C/2 


(3-5) 


Note  that  a  singularity  of  the  integrand  occurs  at  5  =  0,  for  all  values  of  t . 
Paul  (6)  has  presented  a  numerical  technique  for  evaluating  integrals  of  this 
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general  type  taking  into  account  the  singular  point.  The  method  is  essentially 
a  modified  trapezoidal  rule  in  which  weighting  coefficients  are  obtained.  If 
the  range  (0,t)  of  the  integral  is  subdivided  into  subranges  each  of  equal 
length  At,  the  above  integral  for  any  time  t  =  IAt  can  be  represented  as 


<P  = 


(i^l)  jm+l)AT 
m=0  mAT 


fCQdC 

>/  2£ 


(3-6) 


where  f(£)  is  a  polynomial  of  first  degree  in  £  which  approximates  the 
quotient 

J 1+1/2 

as  a  linear  function  over  each  subrange.  The  results  of  the  integration  of 
equation  (3-6)  are  then  put  in  the  form 

(i-D 
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represent  a  set  of  weighting  functions  which  can  be  computed  for  any  interval 
being  considered.  For  large  values  of  m  these  functions  were  written  in  a 
form  more  convenient  for  computational  purposes  by  asymptotic  expansion,  to 
get  the  following  result 
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For  an  illustration  of  the  variation  of  AM  and  EM  functions  with  m,  see 
Fig.  10,  R^  is  termed  the  multiplying  factor  which  for  this  particular 
integral  is  equal  to 


All  of  the  various  integral  expressions  for  the  potential  functions  were 
simplified  to  the  form  of  equation  (5-6)  thereby  permitting  use  of  the  above 
weighting  factors  for  all  integrals.  However,  the  expression  for  the  multi¬ 
plying  factor  varies  for  the  different  integrals. 

The  integration  for  the  outgoing  shear  potential  \|f  is  handled  in  a 
similar  manner,  except  for  a  minor  modification  which  results  from  the  fact 
that  the  shear  wave  travels  outward  with  velocity  c2> 

Consider  for  illustrative  purposes  some  shear  potential  represented 
as 

00 

\|r  =  /  G(t  -  —  cosh  u  )  cosh  uD  du 

J0  c2  d  d  d 

Since  the  incident  wave  strikes  the  shell  at  time  t  =  0,  the  shear  wave  front 
must  also  reach  the  boundary  at  this  instant  in  order  to  satisfy  boundary 
conditions.  Thus  at  t  =  0 

R  cosh  u^  =  R 

and  for  any  time  later 

R  cosh  u^  =  R  +  tc2 

using  the  same  reasoning  as  for  the  outgoing  dilatational  wave.  We  can  express 
the  above  equation  in  non-dimensionalized  form  as 


cosh  u-.  =  1  +  k  t 
2m  c 
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where 

k  .i 

c  cx 

is  a  constant  relating  the  velocities. 

Again  since  is  a  variable  of  integration  we  can  write 

cosh  u  =  1  +  k  £ 

2  c 

where  £  as  in  the  case  of  the  dilatational  wave  is  the  new  variable  of 
integration  vith  limits  0  to  T.  In  terms  of  this  variable  the  integral  is 
now  written  as 

J  0(.  1  ♦  ,  .  5)  VTC  (1  .  kcC)  45 

-  (3-io) 

Note  that  the  shear  wave  starts  from  the  z  axis  at  time  t  =-l/kc  which  is 
earlier  than  for  the  dilatational  wave  since  it  travels  at  a  slower  speed 
(with  fictitious  material  within  .he  shell)  and  must  reach  the  boundary  at 
time  equal  to  zero.  The  numerical  integration  of  equation  (3-10)  is 
accomplished  in  a  similar  method  to  that  which  was  done  for  the  dilatational 
potential.  The  result  can  be  expressed  as 

(i-D 

♦  ■  Z  K  •  S.  •  °1-  *  “«  •  Vi  '  w]  (5-11) 

m=0 


J0 


where  AMm  and  B4m  are  weighting  factors  identical  to  those  derived  previously, 
and  for  this  particular  integral  is 


*  (1  +  W 


1  + 


k  C 
cm 


The  expression  for  is  again  dependent  on  the  form  of  each  integral  considered. 
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For  an  illustration  of  the  variation  of  the  F  and  G  functions  with  time, 
see  Figs.  11  and  12. 


5*3  Solution  of  the  Basic  Equations 

The  equations  of  motion  for  the  shell  are  solved  by  an  iterative 
procedure  known  as  the  Newmark  Beta  Method  (5)  with  which  values  for  the 
potential  functions,  and  the  accelerations,  velocities,  and  displacements 
of  the  shell  are  determined. 

In  general  the  method  consists  of  using  a  step-by-step  integration 
technique  over  successive  time  intervals  assuming  a  specific  variation  of 
acceleration  during  each  interval.  If  we  assume  that  at  time  t  =  (i-l)AT  all 
values  of  potential  functions,  accelerations,  velocities  and  displacements 
are  known,  the  method  becomes  that  of  determining  the  corresponding  values 
at  a  time  t  =  iAt.  For  an  assumed  linear  variation  of  acceleration  over  each 
interval  the  equations  of  interest  are 

a  .  =  a  (a  .  ..  +  a  .) 

n,i  n,i-l  2  n,i-l  n,i' 

At2 

a  .  =  a  .  ,  +  Axa  .  ,  +  -7-  (2H  .  .  +  d  . ) 
n,i  n,i-l  n,i-l  6  '  n,i-l  n,i' 

AT  ••  x  O'12) 

Pn,i  =  &n,i-l  +  2  ^Pn,i-1  +  Pn,i^ 

2 

Pn,i  =  Pn,i-1  +  A^n,i-1  +  ~Z~  ^2^n,i-l  +  ^n,i^ 

where  n  denotes  the  mode  and  i  the  time.  Note  that  in  the  above  equations 

the  velocities  and  displacements  are  expressed  in  terms  of  a  .  .  and  8  ,  , , 

n,i-l  Kn, i-1' 

acceleration  components  which  genexally  are  unknown. 

Values  for  the  acceleration  components  at  time  t  =  IAt  are  assumed 
and  the  velocities  and  displacements  determined  from  the  above  equations;  then 
by  use  of  the  boundary  equations  to  determine  the  values  of  the  potential 
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functions,  values  of  the  acceleration  components  can  he  computed  from  the 
equations  of  motion.  The  computed  values  are  then  compared  with  the  assumed 
values  and  if  they  agree  or  are  within  an  arbitrarily  established  tolerance, 
the  step  is  completed}  if  not,  the  values  of  a n  i  and  0n  ^  just  computed  are 
taken  as  new  assumed  values  and  the  cycle  repeated  until  the  criterion  is 
satisfied.  A  detailed  discussion  of  the  procedure  follows. 

The  equations  of  continuity  at  the  boundary  at  time  t  =  iAr  are 
given  in  equations  (2-49)  and  are  rewritten  here  in  the  form 


LAt 

%,1  ■  F1„<9i>  ♦  f  W  -V 


0 
IAt 


(3-13) 


K,i  -  F2„<91>  ♦  /  «2„(F0'F1— -V  Go'°l--Gi> 


dC 


where  Pln  and  P2n  denote  functions  related  to  the  particle  velocities  in  the 
medium  at  the  boundary  due  to  the  incoming  stress  wave;  their  values  may  be 
determined  directly  at  any  time.  The  integral  expressions  denote  the  con¬ 
tribution  due  to  the  outgoing  waves.  All  values  of  F  and  G  are  known  except 
and  G^^  which  are  to  be  determined  from  the  above  equations.  These  equations 
are  now  written  isolating  the  unknown  values  of  F  and  G  terms 

IAt 

«n,i  *  F1„<81>  *  /  «1„(p0’Fi'  •  •  « 


At 

At 


(3-14) 


+  I0  91n(Fl.l’Fii  « 

IAt 

^n,i  =  P2n^°l^  +  ^  ^n^O'^l'  *  *  *Pi-l'  G0,G1>  ’ '  ,Gi-l^ 

At 

+  f0  Oi-l-V  « 
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The  integrals  with  limits  At  to  IAt  are  evaluated  using  the  numerical 

technique  described  in  a  previous  section. 

The  next  step  is  to  assume  values  for  O  .  and  0  . ,  the  radial  and 

n  ,  i  n  ,  1 

tangential  components  of  acceleration  of  the  shell,  from  which  velocity  and 
displacement  components  are  determined  using  equations  (3-12).  With  these 
values  known,  equations  (3-1*0  become  effectively  two  equations  with  the  two 
unknowns,  F^  and  G^  which  can  then  be  determined. 

The  equations  of  motion  are  given  in  equations  (2-53)  and  are  here 
written  in  the  form 


iAr 

\,l  -S1<°n,i‘  »n,l>  +  PJn<V  +  /  '  •  -ri>  '  V 


dC 


0 
IAt 


(3-15) 


0  .  =  N2(a  p  ,)  +  P4  (O  +  /  Q4  (F  .F. , . . .F. ;  G.G,  ,...G, ) 

Kn,i  x  n,i  Kn,i7  nx  1'  Jn  nN  o7  l7  i7  o7  l7  i' 


from  which  the  acceleration  components  ^  and  i  are  now  computed.  If 
the  computed  values  agree  or  are  within  an  arbitrarily  established  limit  of 
the  assumed  values,  the  step  is  completed;  otherwise  the  computed  values  are 
used  as  the  assumed  values  for  the  next  cycle  of  iteration. 

The  procedure  described  above  requires  that  all  parameters  at  t  =  0 
be  known,  including  the  initial  values  Fq  and  Gq*  Initial  values  of  velocity 
and  displacement  of  the  shell  are  specified  to  be  zero.  From  the  short-time 
approximation  presented  in  Section  3.7  initial  values  of  acceleration  and  G 

o 

were  determined  to  also  be  zero,  and  the  values  of  Fq  to  be  ttie  following 


F 

0 


Ve  . 

IT 


< 


-1 


l  2(-l)n+1 


n  =  0 

(3-16) 

n  =  1,2, . . . 


The  short-time  approximation  also  indicated  that  the  radial  acceleration 
components  near  t  =  0  varied  as  t"^2  and  the  tangential  components  as  t^2; 


thus  to  improve  the  accuracy  of  the  ca] culations,  equations  (3-12)  for  the 
fxrst  step  in  time  were  modified  to 


^3 


2 

3 


a 


'n,l 


a 


n,l 


*n,l  =  5  Pn,l 
Pn,l  =  7  Pn,l 


(3-17) 


Without  this  modification,  the  results  for  the  first  step  cannot  be  brought 
into  acceptable  agreement  with  the  short-time  approximation. 

The  n  =  0  mode  is  simplified  somewhat  because  it  has  no  tangential 
component  of  displacement,  and  contains  only  the  dilatational  component  of 
the  potential  functions. 


3-4  Stresses  of  the  Shell  and  Medium 

Solution  of  the  equations  of  motion  yields  modal  values  of  the 
functions  F  and  G,  and  the  modal  acceleration,  velocity  and  displacement 
components  of  the  shell.  The  stresses  in  the  shell  are  determined  directly 
from  the  displacements  using  equations  (2-55)  and  (2-56). 

Stresses  in  the  medium  at  any  radius  and  time  are  determined  using 
equations  (2-57) •  They  are  here  written  in  the  form 

•0  00  00 

!htr£l  ■  Z  [s„<ei>  *f0  T„<««  *  P  V°>  «]  (3-18) 

p  n=0  u  0 

where  Sn(®1)  represents  the  stress  due  to  the  incoming  wave  and  the  integral 
terms  that  due  to  the  outgoing  waves.  All  are  functions  of  both  radius 
and  time.  Consider  the  stresses  in  the  medium  for  any  radius  equal  to  r  at 
time  t  =  IAt  =  i/N. 

©1  (Fig.  7)  is  now  determined  from  the  equation 
0^  =  arc  cos  ®  (l  -  i/N) 


(3-19) 
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and  takes  on  values  from 


e  =  -  for  t  =  0 
1  r 


to 


91  =  w  for  t  =  (1  + 


The  integral  terms  must  be  recomputed  for  each  radius  and  time  considered. 
As  an  example  consider  the  integral 


<p  = 


-  —  cosh  un )  du, 
ci  1  1 


and  let  i  =  1 —  1,  where  k  =  R/r,  represent  the  nondimens ionalized  time 
k.  r 

r 

after  t  =  0  required  for  the  outgoing  dilatational  wave  front  to  reach  the 

radius  r  being  considered.  Then  for  t  <  the  integral  is  equal  to  zero, 

and  for  t  >  t  the  maximum  value  of  (r  cosh  un )  can  be  written 
r  1 

r  cosh  u.,  =  r  +  t'  c. 


or 


cosh  u,  =  1  +  k  t1 
Im  r 

where  t'  =  T-Tr  represents  the  time  after  the  wave  front  has  reached  the 
radius  r.  We  express  cosh  u^  as 


cosh  u.  ■  1  +  k  t 
1  r* 

where  £  is  the  same  variable  of  integration  as  was  used  previously.  The 
integral  is  then  expressed  in  terms  of  £  as 


<P  = 


F(-1+t  '  -£)  d i 


(3-20) 


The  numerical  method  of  integration  used  is  identical  to  that  described 
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earlier  for  the  integral  at  the  boundary,  with  the  result  that  the  same 
weighting  factors  and  are  applicable.  Note,  however,  that  the  limits 
of  the  integral  are  now  0  to  t';  also  the  multiplying  factor  Rffl  for  this 
particular  integral  becomes 

> IT 


R  = 
m 


1  + 


‘A 


By  a  similar  analysis,  the  integral  of  the  shear  potential 

OP 

♦  =  /  G(t  -  —  cosh  u0)  cosh  Up  du 

J  Cl  co  d  d  d 


can  be  put  in  the  form 


f  -  F  +  T’  '  U  (1  +  k~kJK 

♦  •  /  - 

n 


r  c- 


'J~2t  \jl  + 


kk  J 

r  c 


where 


t*  a  t  -  t 


rs 


For  the  potential  \| r 


c  r 


(3-21) 


The  weighting  factors  AMm  and  EMm  remain  the  same  but  the  multiplying  factor 

for  this  particular  integral  becomes 

V  k  k  (l  +  k  k  &  ) 
re'  r  ctt 


fl  + 


1  *  L 

r  c*m 


46 


3*5  Time  Dependent  Stress  Wave 

The  problem  considered  thus  far  has  dealt  with  an  incoming  stress 
wave  with  a  step  distribution  in  time.  The  results  obtained  for  the  step 
wave  can  be  used  through  the  application  of  Duhamel's  integral  to  find  values 
for  the  response  of  a  shell  and  medium  to  incoming  waves  with  any  time 
variation.  As  an  example,  a  stress  wave  which  decays  exponentially  with  time 
(Fig.  8)  according  to  the  following  equation  is  considered 

-k  ~ 

°P(t)  =  9po(l  ’  e  °  ^5-22) 


represents  the  time  at  which  the  stress  wave  decays  to  zero  and  k  is  a 
parameter  which  is  related  to  the  shape  of  the  curve.  Stresses  at  any  time 
equal  to  IAt  can  be  found  by  the  application  of  Duhamel's  integral,  here 
written  in  the  form 


a(iAr)  =  agt(iAO  + 


iAr 

L 


dt 


ogt(lAr-T)  dT 


(3-23) 


a  t  is  defined  to  be  that  stress  resulting  from  an  incoming  step  wave  of 
amplitude  o^o  in  equation  (3-22).  In  the  numerical  analysis,  the  time 
dependent  wave  is  approximated  by  a  series  of  rectangular  sections  as  illu¬ 
strated  in  Fig.  8. 


3.6  Description  of  the  Computer  Program 

The  problem  was  programmed  modewise  for  a  high  speed  digital 
computer  (CDC  1604)  using  Fortran  language,  and  considering  only  the  first 
three  modes.  Input  data  consist  of  the  following  parameters: 

(l)  The  time  intervals  at  which  computations  are  to  be  performed, 
expressed  as  the  number  of  intervals  required  for  the  incoming  wave  to  travel 
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one  radius  (one-half  transit  time).  The  degree  of  accuracy  achieved  is 
dependent  on  this  parameter,  the  smaller  the  interval  the  more  accurate 
the  results.  However,  the  machine  time  required  for  a  given  number  of  transits 
of  the  incoming  wave  varies  inversely  as  the  square  of  the  interval  size  ap¬ 
proximately,  thus  some  sacrifice  in  accuracy  is  necessary  to  reduce  the  time 
of  computations  required. 

(2)  The  total  time  over  which  the  computations  are  to  be  performed, 
expressed  as  the  total  number  of  time  intervals  to  be  considered.  Generally 
speaking,  all  values  seemed  to  have  reached  their  asymptotic  (static)  values 
within  ten  transit  times  of  the  incoming  wave  across  the  cavity. 

(3)  The  ratio  of  the  moduli  of  elasticity,  E  /E,  where  E  is  the 

s  s 

“plane  strain”  modulus  for  the  shell  and  E  the  modulus  of  elasticity  of  the 
medium. 

(4)  The  mass  ratio,  p  /p,  where  p  is  the  mass  of  the  shell  and  p 

s  s 

the  mass  of  the  medium  per  unit  volume. 

(5)  Poisson's  ratio  of  the  medium. 

(6)  The  ratio  of  the  thickness  of  the  shell  to  its  radius.  When 
considering  a  shell  whose  area,  A,  and  moment  of  inertia,  I,  are  not  directly 
related  to  the  thickness,  A  and  I  must  be  specified  separately. 

(7)  The  amount  of  additional  mass  within  the  shell  expressed  as  a 
fraction  of  the  mass  per  unit  surface  area  of  the  shell. 

(8)  The  number  of  radii  to  which  stresses  in  the  medium  are 

desired. 

(9)  The  time  intervals  at  which  output  data  are  desired. 

(10)  The  angular  increment  at  which  output  data  is  to  be  computed. 


Because  of  symmetry  only  values  between  0  and  180  degrees  need  be  considered. 
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Additional  input  quantities  in  the  case  of  the  exponentially- 
decaying  stress  wave  are  tq  and  k,  parameters  which  define  the  duration  and 
shape  of  the  stress  pulse,  respectively. 

Output  data  consist  of  the  acceleration,  velocity,  displacement, 
and  stress  components  of  the  shell  for  specified  angles  and  times;  and  the 
stresses  in  the  medium  for  specified  radii,  angles  and  times. 

3.7  Short-Time  Approximation 

As  a  check  on  the  accuracy  of  the  machine  solution  for  short  times, 
the  boundary  equations  and  the  equations  of  motion  were  solved  approximately 
by  making  a  series  expansion  of  all  pertinent  functions  in  terms  of  time  as 
the  independent  variable.  Although  the  following  discussion  is  limited  to  the 
case  of  the  incoming  dilatational  wave,  the  basic  principles  are  the  same  for 
either  of  the  types  of  wave  considered. 

The  expressions  for  the  velocities  and  stresses  in  the  medium 
around  the  boundary  due  to  the  incident  wave  all  have  been  written  thus  far 
in  terms  of  0^,  the  position  angle  of  the  wave.  In  terms  of  nondimensionalized 
time  t,  all  functions  of  9^  can  be  written  in  terms  of  x  using  the  following 
relations 


cos  0^  =  1  -  x 


The  integral  values  representing  the  effects  of  the  outgoing  waves 
are  also  expressed  as  functions  of  x.  The  following  example  will  illustrate 
the  technique  used  to  accomplish  this.  For  example,  consider  the  transforma¬ 


tion  of  the  integral 


^9 


qp  =  /  P(t  -  —  Cosh  u. )  du.. 

J  0  C1  1  x 


As  was  shown  in  Section  3-2,  cosh  in  the  integral  varies  in  value  for  any- 
given  time  t  from  1  to  1  +  t.  It  is  convenient  for  purposes  of  analyzing  the 
integral  to  represent  the  variable  cosh  u^  as  follows 

cosh  u^  »  1  +  wt  (3-25) 

where  t  is  now  a  fixed  value  in  the  integration  and  u  is  defined  as  the 
variable  of  integration  whose  value  ranges  from  0  to  1.  The  function 

F(t  -  —  cosh  u.. )  =  F(-l  +  t(1-u>)) 

C1 

is  expanded  in  terms  of  a  power  series  as 

00 

F(-l  +  t(l-u>) )  =  Fo  +  7iT1(l-u))i  (3-26) 

i=l 

where  7 ^  are  unknown  coefficients  of  the  series.  Since 


the  integral  can  now  be  written  as 


from  which  after  performing  the  indicated  integration 

<p  [Fo  +  T(j  71  -  %  F0)  +  t2  (•••  ]  (3-27) 

The  integrals  which  represent  the  effects  of  the  outgoing  shear 
wave  can  also  be  transformed  in  somewhat  similar  manner.  Consider  for  example 
the  integral 
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OS 

=  /  G(t  -  cosh  uj  cosh  u, 

J  n  co  c- 


2  du2 


In  this  case  cosh  varies  from  1  to  1  +  k  \  where  k  =  —  and  t 
c  c  c  c^ 

Therefore  we  write 


tc. 


cosh  Ug  =  1  +  cokj  (3-28) 

where  w  is  now  the  variable  of  integration  ranging  in  value  from  0  to  1  in 
the  integral.  The  function 


G(t  -  £-  cosh  u?)  =  G(-  ~-  +  t(1-u))) 
C2  *  c 


is  expanded  in  terms  of  a  power  series  as 

00 

0(“  +  t(1-u))  =  Go  +  £  eiTi(l-u)i 


i=l 


where  €1  are  unknown  coefficients  of  the  series.  From  equation  (3-28) 


>/  2k  t 


du„  = 


2 

V  CO 


t  r  uk  i 

=  1  -  T"  ♦  #  <“V>  -  ••• 

CO  * 


dw 


The  integral  written  in  terms  of  the  variable  w  is  now 


'/k  t 


/0  ^[°o  +  V(1’w)+***]  +  *T“  +  •••]  l1  +  «v]  dw 


■f2  '■'0  >/w 

and  performing  the  integration 


*  ■  ^7  [t(f  ao>  *  6l  +  B  kc  =o>  +  '5<  •  •  •  ]  (3-29) 


The  same  basic  technique  is  applied  to  all  the  integrals  so  that  now  the 
effects  of  both  the  incoming  and  outgoing  waves  can  be  written  in  terms  of  t. 
Substitution  into  the  continuity  equations  and  the  equations  of  motion  yield 
the  following  for  any  mode 
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l  : 


a  = 

^  =  £  f2i(F,G)t  2 
i=l 

(3-30) 

a  +  =  £  f31(F,G)x  2 

1=1 

*  i  _  1 

P  +  t2±(a,p)  =  £  2 

i=l 

where  fl^,  f2^,  f3^,  and  f4^  are  coefficients  which  contain  certain  elements 

F  ,  G  ,  7.  and  e.  of  the  potential  functions;  21  and  22  are  known  functions 
0  0x1 

of  the  displacements. 

To  solve  the  above  equations,  the  displacement  components  of  the 
shell  are  expressed  as  Frobenius  (8)  type  series 


fl^Gb 


to 


i=l 


(3-31) 


where  p^,  q^>  and  c  are  unknown  coefficients.  Substitution  into  equations 
(3-30)  results  in  the  following  set  of  equations 

00  *  l 
^  Pj_  •  (c+i)  •  Tc+i_1  =  £  fl^Gb  2 
i=l  i=l 


00 

I 

1=1 


(c+i) 


c+i-1 


^  f2i(F,G) 


i 

T 


i=l 


1 

2 


(5-32) 
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OD 

v  p1  *  (c+i)-(c+i-i) 


,0*1-2  + 


i=l 


£  ll(pi,qi)'iC+i  =  £  f31(F,G)x 
^  i=l 


i=l 


00  00  op 

£  •  (c+i)-(c+i-l)  •  Tc+i'2  +  £  /2(pi>qi)TC+i  =  £ 

i-1  i=l  i=l 


The  coefficient  c  is  now  determined  by  inspection.  Then  through  a  step  by 
step  process  which  involves  the  equating  of  coefficients  of  like  powers  of  t , 
values  of  p^  q^,  7 and  e i  are  determined.  The  following  equations  are  then 
used  to  find  values  of  the  potential  functions,  and  the  displacement  components 
of  the  shell. 


f(t)  =  fo  +  ^ 
i=l 
00 

C(t)  -  Go  +  £  €1T1 
i=l 

»  (3-33) 

<*(■0  =  Yj  Pi  TC+i 

i=l 

00 

P(l)  =  ^  1iTC+i 
i=l 


Velocity  and  acceleration  components  may  be  determined  by  differentiation  of 
the  above  expressions  for  the  displacements. 
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CHAPTER  IV 

*.  DISCUSSION  OF  RESULTS 


4.1  General 

Results  of  computations  performed  to  determine  the  effect  of  the 
various  parameters  are  discussed  in  this  chapter. 

Although  equations  presented  throughout  the  study  have  been  written 
to  include  an  infinite  number  of  modes,  the  greater  part  of  the  actual  calcu¬ 
lations  performed  and  presented  here  are  the  results  obtained  considering  only 
the  modes  n  =  0,  1,  and  2.  It  is  important  to  note  that  during  envelopment  of 
the  shell  by  the  plane  stress  wave,  a  Fourier  series  representation  of  the 
incoming  wave  is  objectionable  in  that  the  series  at  this  stage  is  slowly  con¬ 
vergent,  thus  necessitating  a  large  number  of  modes  to  accurately  represent  the 
plane  wave.  However,  after  passage  of  the  wave  across  the  cavity,  the  Fourier 
expansion  of  the  incoming  stresses  around  the  boundary  results  in  coefficients 
of  all  modes  except  n  =  0  and  2  becoming  identically  equal  to  zero  for  the 
plane  dilatational  wave,  and  coefficients  of  all  modes  except  n  =  2  becoming 
identically  equal  to  zero  for  the  plane  shear  wave.  Therefore,  stresses  due 
to  the  outgoing  waves  in  modes  corresponding  to  those  of  the  incoming  wave 
whose  coefficients  become  zero  must  also  eventually  vanish  at  long  times.  The 
limited  study  conducted  for  modes  greater  than  n  =  2  indicated  that  the 
maximum  effect  of  the  higher  modes  occurs  within  one  transit  time  of  the 
incident  wave  and  rapidly  decays,  thus  contributing  relatively  little  to  the 
maximum  response  of  the  shell  which  occurs  after  several  transit  times. 

However,  for  determining  the  early  time  response  of  the  shell  and  medium, 
the  higher  modes  are  significant  and  should  be  considered  in  further  extension 


of  this  work. 
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In  the  tables  and  figures  to  be  discussed,  quantities  given  in 

non-dimensionalized  units  are  defined  by  equations  (2-36)  and  (2-50).  Stresses 

are  given  in  units  of  the  absolute  value  of  the  amplitude  (|o  |  or  |a  |)  of  the 

P  S 

incident  wave;  a  negative  stress  means  a  compressive  response  to  an  incoming 
compressional  (Fig.  l)  or  a  positive  shear  wave  (Fig.  2).  The  physical 
properties  of  the  shell  relative  to  those  of  the  medium,  as  well  as  the  nature 
of  the  incoming  wave  are  indicated  on  the  graphs.  Unless  otherwise  stated, 
the  shell  is  considered  to  he  an  unstiffened  one  so  that  its  cross  sectional 
area  and  moment  of  inertia  are  related  to  the  thickness  as  given  hy  equations 
(2-52) •  Except  where  indicated,  there  is  assumed  to  he  no  additional  mass 
within  the  shell.  Numeral  subscripts  denote  the  mode  number. 

The  shell  and  medium  have  been  assumed  to  exhibit  linearly  elastic 
behavior  throughout  their  stress  histories,  which  for  the  practical  problem 
does  not  permit  evaluation  of  any  spalling  or  non-elastic  effects. 

Values  of  stresses  given  are  in  addition  to  those  which  exist  prior 
to  the  arrival  of  the  incident  wave.  For  the  elastic  case,  the  effect  of 
prior  stresses  such  as  those  resulting  from  the  overburden  may  be  taken  into 
account  by  merely  adding  them  to  stresses  caused  by  the  incident  wave. 

For  clarity  of  presentation  and  because  of  the  impracticability  of 
including  solutions  for  all  possible  permutations  of  the  parameters  involved, 
the  discussion  in  this  chapter  is  limited  to  a  few  representative  cases. 

4.2  Modal  Response  of  the  Shell  and  Medium 

Figures  11  and  12  illustrate  the  shape  of  the  modal  components  of 
the  dilatational  and  shear  potentials  obtained  in  the  solution  to  a  typical 
problem.  It  appears  that  a  singularity  occurs  at  one  transit  time  in  the 
case  of  the  F  functions  resulting  in  the  s Light  irregularity  of  the  curves  at 
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this  point.  Since  computed  iralues  of  stresses  were  determined  to  be  rather 
insensitive  to  relatively  large  variations  in  values  of  the  F  and  0  functions, 
the  effect  of  the  irregularity  would  seem  to  be  slight. 

Figures  13  through  20  show  modal  acceleration,  velocity,  displace¬ 
ment  and  stress  components  for  the  shell  and  stress  components  for  the  medium 
at  the  boundary,  as  they  vary  with  time.  Static  values  shown  were  computed 
using  the  method  given  in  Appendix  B. 

The  high  accelerations  computed  near  the  beginning  are  not  truly 
representative  of  the  actual  case,  since  they  are  the  result  of  assumptions 
made  earlier  in  deriving  the  equations  for  the  shell.  The  shell  was  repre¬ 
sented  by  a  line  describing  its  middle  surface  which  permits  no  variation  in 
accelerations,  velocities,  and  displacements  of  particles  through  the  actual 
thickness.  Also,  no  provision  was  made  for  refraction  of  the  incident  wave 
through  the  shell  lining.  These  limitations  restrict  the  applicability  of 
the  solutions  to  a  shell  whose  thickness  is  small  relative  to  its  radius. 

The  n  =  1  mode  is  primarily  a  translational  one  which  accounts  for 
the  rigid  body  translation  of  the  shell  after  it  has  been  enveloped  by  the 
incident  wave.  Thus,  it  can  be  seen  that  the  velocity  components  for  this 
mode  approach  constant  values  equal  to  the  velocity  of  the  medium  behind  the 
incident  wave  front,  and  displacements  grow  without  bound  reaching  a  straight 
line  variation  with  time.  Note  that  the  s^resses  contributed  by  this  mode 
reach  their  peak  values  within  one  transit  time  and  quickly  damp  out, 
approaching  zero  asymptotically.  For  the  incident  shear  wave,  the  n  =  0  mode 
is  also  a  rigid  body  movement  which  accounts  for  rigid  body  rotation,  and 
contributes  little  to  the  stresses. 

Modal  quantities  obtained  are  coefficients  of  Fourier  series; 
therefore,  the  total  response  or  effect  is  determined  by  adding  the  coefficients 
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multiplied  by  the  appropriate  sine  and  cosine  terms  for  any  desired  angle. 
Figures  21  through  24  show  the  time  variation  of  stresses  in  the  shell  and 
medium  for  various  angles  when  the  first  three  modes  are  summed.  The  maxi¬ 
mum  stress  in  the  shell  for  any  angle  may  be  determined  by  adding  the  bending 
stress  to  the  hoop  stress.  This  is  indicated  in  Figs.  21  and  24  by  the 
dotted  line  above  the  hoop  stress. 

Figure  25  is  given  to  illustrate  the  relative  magnitudes  of  the 
hoop  stresses  in  the  shell  and  medium  for  several  thicknesses  of  shell.  This 
also  shows  the  effect  of  varying  the  relative  thickness  of  the  shell  on  the 
hoop  stress  in  the  medium.  The  dotted  line  indicating  the  hoop  stress  in 
the  medium  for  an  unlined  cavity  was  obtained  from  the  report  by  Paul  (6). 

Figures  26  and  27  show  how  the  relative  thickness  of  the  shell 
affects  the  radial  and  shear  stresses  in  the  medium. 

As  was  discussed  earlier,  stresses  in  the  medium  for  any  radius 
can  be  determined  by  reevaluating  the  integral  terms  which  represent  the 
effects  of  the  outgoing  waves,  and  adding  them  to  the  Fourier  expansion  of 
the  incident  wave.  Figure  28  shows  the  modal  and  total  radial,  hoop,  and 
shear  stresses  which  were  computed  for  a  time  equal  to  10  transit  times. 

These  values  are  compared  later  with  the  static  stresses,  but  on  this  figure 
the  static  stresses  do  not  differ  by  more  than  the  thickness  of  the  lines, 
and  therefore  are  not  shown.  The  time  variation  of  the  radial  and  hoop 
stresses  in  the  medium  for  various  radii  are  shown  in  Fig.  29  for  the  incident 
dilatational  wave,  and  in  Fig.  30  for  the  incident  shear  wave. 

4.3  Short-Time  and  Asymptotic  Comparisons 

A  method  for  obtaining  a  solution  to  the  problem  which  is  ac^^ate 
for  very  short  times  (t  «  l)  was  presented  in  Section  3*7*  This  was  desirable 
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to  validate  the  machine  solution  and  to  determine  the  effect  of  the  size  of 
time  interval  selected.  The  results  for  a  representative  problem  are  shown 
in  Figs.  31  through  33  for  the  incident  dilatational  wave,  and  in  Figs.  3^ 
through  36  for  the  incident  shear  wave.  Four  terms  of  the  series  representing 
the  F  and  G  functions,  and  three  terms  for  other  quant it ies  were  used  in  the 
short-time  solution. 

As  can  he  seen  from  the  graphs,  good  agreement  was  obtained  for 
very  small  values  of  time,  somewhat  shorter  time  being  obtained  for  the  shear 
wave  as  compared  to  the  dilatational  wave.  The  shorter  time  results  from  the 
nature  of  the  forcing  functions  (Eqs.  2-19  and  2-32)  which  indicate  a  more 
rapid  rise  in  the  incoming  stresses  for  the  incident  shear  wave. 

Within  the  range  of  time  for  which  the  short-time  solution  is  valid, 
decreasing  the  size  of  time  interval  for  each  step  of  the  machine  solution 
resulted  in  closer  agreement  between  the  two  methods,  as  is  to  be  expected. 

It  also  indicated  that  the  stresses  and  displacements  are  not  as  sensitive  to 
variations  in  the  interval  size  as  are  the  F  and  G  functions. 

At  the  other  end  of  the  time  scale,  i.e.,  at  a  relatively  long  time 
after  passage  of  the  incident  wave  front  across  the  cavity,  another  check  on 
the  accuracy  of  the  machine  solution  is  afforded  by  the  asymptotic  approach 
of  all  values  to  the  static  results.  Figures  presented  thus  far  have  shown 
that  the  static  condition  is  approached  well  within  ten  transit  times. 

The  time  interval  used  in  the  machine  solution  affects  the  stability 
of  the  results  for  long  times.  This  is  indicated  in  Fig.  37  which  shows  the 
variations  in  computed  values  of  the  displacement  components  for  mode  n  =  2 
at  relatively  long  times,  for  different  time  intervals.  N  represents  the 
number  of  time  steps  required  for  the  incident  wave  to  travel  a  distance  equal 
to  the  radius  of  the  opening,  and  l/N  defines  the  interval  size.  Note  the 
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smaller  graph  which  shows  the  percentage  difference  between  the  computed 
values  at  ten  transit  times  and  the  static  values. 

As  the  interval  is  decreased  the  machine  solution  at  long  time 
approaches  the  static  solution  more  closely.  Below  a  certain  size  of  time 
interval,  there  is  little  difference  in  the  results,  which  indicates  asymp¬ 
totic  convergence  to  the  correct  solution.  For  this  particular  problem, 

N  -  30  seems  to  be  "critical"  in  that  for  N  <  }0,  wide  variations  in  computed 
values  occur.  As  the  mode  number  increases,  the  "critical"  value  of  N  in¬ 
creases  rapidly,  and  the  requirements  of  computer  storage  and  calculation 
time  become  decisive  factors  which  make  impractical  the  study  for  long  times 
of  modes  much  larger  than  2.  The  interval  size  selected  for  all  problems 
solved  (exclusive  of  the  study  to  determine  the  effect  of  the  interval  size) 
was  set  equal  to  l/40  (N  =  40)  of  the  half  transit  time  of  the  incident  wave. 

For  the  static  case,  only  the  modes  n  =  0  and  2  yield  values  other 
than  zero.  Table  1  compares  values  of  shell  stresses  and  displacement  com¬ 
ponents  in  these  modes  obtained  from  the  computer  solution  to  a  particular 
problem  at  a  time  equal  to  10  transit  times,  with  the  static  solution.  Most 
pairs  of  values  differ  by  less  than  one  percent.  Comparable  agreement  of 
stresses  in  the  medium  at  various  radii  are  shown  in  Table  2.  Values  in 
Tables  1  and  2  were  obtained  from  the  solution  to  a  problem  whose  parameters 
were:  tje  =  1+.0,  np  =  5*0,  rjt  =  .05,  and  v  =  .25. 

Although  the  figures  and  tables  presented  above  were  for  a  particular 
problem,  the  discussion  given  is  applicable  to  all  problems  which  were  solved. 
Changing  the  physical  characteristics  of  the  shell  and  medium  within  the  range 
of  values  studied  had  hardly  any  effect  on  the  degree  to  which  the  long  time 
machine  solution  and  the  static  solution  agreed.  Also,  the  absolute  value  of 
all  quantities  which  should  asymptotically  approach  zero  became  less  t.ian 
.00005  in  each  case  well  before  ten  transit  times. 
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4.4  Effect  of  Parameters 

Studies  were  conducted  to  determine  the  effect  of  each  of  the 
following  parameters : 


r\t  =  -  V  *  V  (medium) 

for  an  unstiffened  shell  without  additional  mass.  Since  it  was  impractical 
to  take  into  account  all  permutations  of  the  above  parameters,  a  basic  shell 
where 

\  =  4.o  np  -  3.o 

\  =  -05  v  =  .25 

was  considered  from  which  each  parameter  was  separately  varied  to  determine 

its  effect  on  the  resulting  stresses  and  displacements.  Calculations 

performed  were  only  for  the  case  of  the  incident  dilatational  (P)  wave. 

Of  particular  interest  was  the  determination  of  the  maximum  dynamic 

stresses  and  displacements  (not  including  rigid  booty  translation)  due  to  the 

incident  stress  wave.  Tables  3  through  6  compare  the  maximum  values  obtained 

in  the  machine  solution  with  the  static  solution.  "DLF,"  termed  the  dynamic 

load  factor,  is  defined  as  the  factor  by  which  the  displacement  or  stress 

produced  by  a  applied  as  a  static  load  should  be  multiplied  in  order  to  obtain 
P 

the  maximum  dynamic  value.  Figures  38  through  41  are  given  to  graphically 
illustrate  the  variations  in  stresses  and  displacements,  both  static  and 
dynamic,  in  the  range  of  parameters  considered.  Stresses  are  given  in  units 
of  | 0^ |  and  displacements  in  non-dimens ionalized  units  defined  by  equations 
(2-36)  and  (2-50). 
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Table  3  and  Fig.  38  illustrate  the  effect  of  increasing  the 
relative  thickness  of  the  shell.  The  range  of  values  selected  for  r\  is 

T> 

probably  much  greater  than  is  practical  or  justified  by  the  assumptions  of 
the  analysis,  and  was  considered  only  to  determine  the  trend  of  the  results. 
As  t)^  increased,  displacements  and  hoop  stresses  of  both  the  shell  and 
medium  decrease;  however,  the  bending  stress  of  the  shell  and  the  radial  and 
shear  stresses  in  the  medium  at  the  boundary  increase.  Figure  25  shows  the 
time  variation  of  the  hoop  stress  in  the  shell  and  medium  for  various  thick¬ 
nesses  of  shell,  including  the  case  of  the  unlined  cavity. 

Increasing  results  in  a  rapid  increase  in  the  displacements  and 
stresses  of  the  shell  and  a  much  lower  rate  of  attenuation  of  the  hoop  stress 
in  the  medium.  See  Fig.  39  and  Table  4. 

Table  5  and  Fig.  40  show  that  the  parameter  has  no  effect  on  the 
static  results  but  does  affect  the  dynamic  response.  Increasing  r\  increased 
the  maximum  response  in  both  the  medium  and  shell. 

Figure  41  and  Table  6  show  that  within  the  range  of  Poisson's  ratio 
for  the  medium  considered,  as  v  increases  displacements  and  stresses  in  the 
shell  decrease,  with  very  little  additional  reduction  of  the  medium  hoop 
stress . 


Additional  mass  within  the  shell  was  shown  to  contribute  relatively 
little  to  the  overall  response.  This  is  mainly  due  to  the  simplifying  assump¬ 
tions  which  were  made  in  deriving  the  equations  of  motion  to  account  for  the 
additional  mass.  Figure  42  shows  the  effect  of  the  presence  of  additional 
mass  equal  to  40tt  times  the  mass  of  the  shell  on  the  dynamic  response  of  a 
particular  shell.  Additional  mass  decreases  the  displacements  and  increases 
the  dynamic  stresses  in  mode  n  =  1. 
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4.5  Response  to  Time  Varying  Incident  Wave 

The  results  obtained  from  a  solution  to  an  incident  wave  with  a  step 
variation  in  time  was  shown  to  be  useful  by  application  of  Duhamel's  integral  in 
determining  the  response  to  any  time  varying  stress  wave.  The  case  of  the 
exponentially  decaying  wave  was  considered,  the  results  of  which  are 
illustrated  in  Fig.  45 .  Substantial  reductions  in  the  maximum  stresses  can 
be  expected  as  the  duration  becomes  smaller.  In  Fig.  44  is  shown  the  effect 
of  a  linear  rise  in  the  amplitude  of  the  incident  wave.  Note  that  for  a 
wave  with  a  linear  rise  followed  by  a  step  variation  in  time,  veiy  little 
decrease  in  maximum  stresses  occurred. 

4.6  Comparison  with  Previous  Work 

Baron  (l),  using  a  different  method  of  analysis,  investigated  the 
dynamic  response  of  two  shells  with  different  physical  characteristics  sub¬ 
jected  to  an  incident  plane  dilatational  wave.  Figures  45  through  49  compare 
his  results  for  the  modal  values  of  stresses  and  displacements  of  the  thin 
shell,  with  results  obtained  by  the  method  of  solution  given  in  this  study. 
Similar  results  were  obtained  considering  his  so-called  stiff  shell.  Although 
the  shape  of  the  response  curves  can  be  considered  similar  in  the  two  reports, 
the  magnitude  of  the  dynamic  response  in  Ref.  (l)  seems  consistently  higher 
than  in  the  present  report,  and  the  long  time  results  asymptotically  approach 
values  higher  than  the  static  solution.  The  following  table  compares  the 
maximum  modal  stresses,  for  Poisson's  ratio  of  0.25> obtained  in  the  two  reports, 
and  also  shows  values  for  the  static  case.  Stresses  are  given  in  units  of 
|a  |.  The  bending  stresses  for  the  stiffened  shell  cannot  be  compared  directly 
since  d,  the  distance  from  the  neutral  axis  of  the  shell  to  its  extreme  fiber, 
is  not  stated  in  Ref.  (l). 
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Quant. 

Mode 

Thin  Shell 

Stiff  Shell 

Maximum 

Stresses 

Static 

Maximum 

Stresses 

Static 

Ref. 

(1) 

Current 

Work 

Ref. 

(1) 

Current 

Work 

a99 

0 

-4.79 

-4.56 

-4.10 

-4.79 

-4.08 

-3.67 

a99 

2 

4.69 

4.11 

3.67 

3-42 

3.20 

2.95 

°sb 

0 

0.043 

0.043 

0.039 

0.55 

4.08d 

3-673. 

asb 

2 

0.140 

0.137 

0.122 

1.70 

13-29d 

11.89d 

4-7  Conclusions 

Conclusions  drawn  from  the  results  of  the  analysis  are: 

(1)  The  method  which  has  teen  presented  is  practical  for  effective 
computation  of  the  dynamic  response  of  a  cylindrical  shell  embedded  in  an 
elastic  medium  when  subjected  to  plane  dilatational  or  shear  waves .  The 
solution  presented  herein  is  believed  correct  since  it  was  checked  by  independent 
methods  at  short  and  long  times. 

(2)  Peak  stresses  and  displacements  in  both  the  medium  and  shell 
occur  sometime  after  the  transit  of  the  incident  wave  across  the  cavity;  of 
the  problems  solved,  both  the  average  and  mean  time  at  which  the  peak  values 
occurred  was  equal  to  3  transit  times.  The  dynamic  effect  measured  in  terms 
of  the  ratio  of  the  maximum  stresses  and  displacements  to  the  static  values 
varied  within  a  relatively  small  range.  The  average  value  of  this  ratio  was 
1.12  and  the  mean,  1.11.  Thus,  for  the  cases  considered  at  least,  the  maximum 
stresses  and  displacements  to  be  expected  for  a  particular  situation  can  be 
roughly  approximated  by  determining  the  static  values  and  by  multiplying  them 
by  a  factor  of  1.1. 
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(3)  The  largest  stresses  in  the  shell  occurred  for  a  relatively 
thin  liner  in  a  medium  with  a  low  modulus  of  elasticity,  low  mass  density, 
and  low  Poisson* s  ratio.  The  greatest  reduction  of  the  hoop  stress  in  the 
medium  as  compared  to  the  unlined  cavity  results  from  a  relatively  thick 
liner  in  a  medium  with  low  modulus  of  elasticity  and  high  mass  density. 
Additional  mass  within  the  shell  has  relatively  small  effect  on  the  dynamic 
stresses. 

(4)  The  practical  value  of  tunnel  linings  to  reduce  the  maximum 
stresses  in  the  medium  depends  on  the  several  conditions  mentioned  above,  and 
on  the  magnitude  of  the  incident  stress  wave.  Under  favorable  conditions, 
reduction  of  stresses  on  the  order  of  30  percent  or  more  is  possible.  However, 
for  materials  such  as  granite,  smaller  reductions  can  be  expected  for  steel 
liners  of  practical  size. 

(5)  Certain  assumptions  made  concerning  the  behavior  of  the  shell 
have  limited  th  applicability  of  the  analysis  to  relatively  thin  liners. 

Future  studies  of  the  behavior  of  thick  shells  would  indicate  the  effect  of 
the  approximations  used  herein.  The  analysis  has  also  been  based  on  the 
assumption  that  the  behavior  of  both  the  liner  and  medium  is  linearly  elastic 
at  all  times.  Perhaps  a  more  desirable  condition  would  be  one  in  which  some 
inelastic  behavior  is  permitted  to  take  place  in  the  medium  surrounding  the 
shell  or  in  the  shell,  or  one  in  which  some  inelastic  energy  absorbing  medium 
such  as  cinders  or  foamed  plastic  surrounds  the  shell.  As  a  subject  of  future 
study,  it  is  recommended  that  the  behavior  of  thick  shells,  and  thin  shells 
surrounded  by  some  energy  absorbing  layer,  be  considered. 
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APPENDIX  A 

DERIVATION  OF  THE  EXPRESSIONS  INVOLVING  THE  POTENTIAL  FUNCTIONS 


A.l  General  Form  of  the  Potential  Function 

The  dilatational  potential  function  must  satisfy  the  wave 

equation 


isL+  1 

dr2  r 


i  a2cp  i  a2<p 


2  ^2 
r  o9 


2  *.2 
c1  ot 


(A-l) 


9  in  the  case  of  the  incoming  dilatational  wave  is  an  even  function  of  Q  and 
thus  can  be  expressed  as  a  cosine  Fourier  series 


9  =  fn^r't^  cos  nQ  (A-2) 

n=0 

where  f  is  a  function  of  r  and  t  representing  the  modal  coefficient  of  the 

potential  function.  By  substitution  of  the  above  expression  into  equation 

(A-l)  we  obtain  the  equation  that  must  be  satisfied  by  f  . 

n 


d2f  .  Sf  2  _  d2f 

_ _ n  1  n  n  r  _1 _ n 

x  2  r  ST  “  2  n  =  2  ^2 
or  r  c^  dt 


(A-3) 


The  general  solution  of  f  is  assumed  to  be  of  the  form 

n 

f  »  rn  R  (r,t) 
n  nx  9  ' 

Then  Rn  is  some  function  of  r  and  t  satisfying 


d\ 


2^ 


(2n-H) 

r 


(A-4) 


(A-5) 


If  Rn  is  a  solution  it  can  be  shown  that  the  corresponding  equation  for 
Rn+1  is  satisfied  by 


n+1 


,  dR 

1  n 

r  Sr” 


(A-6) 


By  repeated  application  of  this  result  it  can  be  shown  that  equation  (A-5) 
is  satisfied  by 
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R  =  (-  I-)  R 
n  vr  dr'  « 


where  R  is  the  solution  of 
o 


(A-7) 


d2R  0  rd2R  ,  dR  -i 

o  2  |  o.l  o  I 


(A-8) 


The  solution  to  this  equation  for  the  case  of  a  wave  diverging  from  a  center 
is 


00 

R„  =  /  F(t  -  —  cosh  un )  du., 

°  J0  c!  1  1 


(A-9) 


For  the  proof  of  this,  see  Lamb  (4).  Therefore,  the  coefficient  f  of  the 

n 

potential  function  is  written  as 


op 

fn  =  fQ  F(t  -  fj;  cosh  V  dul 


For  modes  0,  1,  and  2,  this  can  be  written 


(A-10) 


00 

f»  ■  L 


cosh  u^  du^ 


cosh^u-, 


du^  + 


cosh  u 


du 


1 


where  primes  indicate  the  derivative  of  the  function  F  with  respect  to  its 

argument.  Note  that  the  expression  for  the  n  =  2  mode  contains  derivatives 

of  two  different  orders.  It  is  convenient  to  express  the  function  f  in 

n 

terms  of  derivatives  of  a  single  order  which  may  be  done  through  integration 


6  7 


by  parts.  For  the  n  =  2  mode  we  then  wish  to  change  the  second  integral  to 
an  integral  involving  the  second  derivative  of  the  function  F.  Integrating 
by  parts  gives 


u,  u, 

n lm  n  1m 

/  =  £*t) 

J  0 


J  o 


p\Lm 

/  n-ds 
J  0 


where 


Thus 


I  =  5" (6X) 

d6  =  -  ~  F"(61)  sinh  ^  di^ 


Tj  =  sinh  u. 


dr)  =  cosh  u^  du^ 


00  00 

f0  F'<51>  c°sh  “1  d"l  -  IT  f0  F'<h> 


sinh  u^du^  + 


F,(61)  sinh  ux 


lm 


The  second  term  can  be  shown  to  equal  zero  since  for  its  lower  limit 
sinh  u1  «  0;  and  for  its  upper  limit,  F'  is  the  integral  of  the  function  F11 
at  the  wave  front  which  under  the  assumed  initial  conditions,  does  not  exist. 
Therefore,  the  coefficient  of  the  potential  function  for  the  n  =  2  mode  can 
be  written  as 


00 

=  —r  /  F"(t  -  —  cosh  u, )  cosh  2u,  du. 
cf  J  0  C1  1  J 


(A-ll) 


By  a  similar  process  the  coefficients  of  the  potential  function  for  any  other 
mode  can  be  reduced  to  the  form 


/  >n  r» 

f  =  *—  *  /  F°(t  -  —  cosh  u, )  cosh  nu,  du. 

n  „n  JA  c.  1  11 


(A-12) 


The  shear  potential  function  must  satisfy  the  wave  equation 

djjr  1  by  1_  dj|r  _  _1_  dft 
Sr2  +  r  57  +  r2  *  c2  dt2 


(A -13) 
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^  in  the  case  of  the  incoming  dilatational  wave  must  be  an  odd  function  of  9 
and  can  be  represented  as 

00 

♦  =  ^  gn(r,t)  sin  n0 
n=l 

and  proceeding  exactly  as  in  the  case  of  the  dilatational  potential  the 
general  expression  for  the  coefficient  g^  is  found  to  be 

00 

g  =  /  Gn(t  -  cosh  U0)  cosh  nup  du  (A-14) 

n  cl  J0  C2  d  d 


A. 2  Velocity  Terms 

The  equations  for  the  velocity  components  of  a  particle  on  the 
boundary  due  to  the  outgoing  waves  are 


(A -15) 


The  components  of  velocity  and  the  potential  functions  are  expanded  in  series 


u  =  X  *„<r’t)  cos  n0 

n=0 

00 

v  =  £  bn(r,t)  sin  n0 
n=l 

oo  (A -16) 

‘P  =  Yj  fn(r,t)  cos  n0 
n=0 


00 


n=0 


sin  nQ 
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Substitution  into  equations  (A-15)  gives 

.  ,  .  .  ^fn  n  dgn 

an(out)  =  S5t  +  FST 

df  b2g 
bn(out)  =  -  7  ‘  55t 


(A-17) 


The  expressions  for  the  coefficients  f^  and  gn  are  given  in  equations  (A-12) 
and  (A-14).  Using  these  in  the  above  equations  and  by  the  application  of 
integration  by  parts  we  get  the  following  velocity  terms 

00 

an(out)  =  -  J  cosh  u^  cosh  nu^  du^ 

c 1  ° 

00 

+  L  Gn+2(42}  sinh  u2 sinh  nu2  du2 


00 

bn(out)  =  -  111 ■  J  Fn+^(|^)  sinh  u^  sinh  nu-^  du^ 


(A-18) 


00 

+  f 

n+1  J  a 


Gn+^(6g)  cosh  u2  cosh  nu^  du^ 


A, 3  Stress  Terms 

Stresses  in  the  medium  due  to  the  outgoing  waves  are  written  in 
series  form  as 


rr 


00 

*  £ 


t)  cos  n& 


n=0 


a  =  )  B  (r,t)  sin  n£ 

TQ  Lx  Vl  * 

n=l 

00 

°ee  =  1  cn(r’t>  cos  n0 
n=0 


(A-19) 
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By  substitution  into  equations  (2-LOj  the  coefficients  A  ,  B  ,  and  C  are 
J  n  n  n 


expressed  in  terms  of  the  coefficients  of  the  potential  functions,  f^  and  g^ 

2 

rb  f  dg  ng  i 

An(out)  -  (X  ♦  *)  *  a  5?  -  -§] 

v  rSf  2 

.  K  I  n  ,  n  n  „  &n 

+  7  |5T  +  7  gn  ‘  T  fn  ‘  n  ST] 

I,  (  *  ^  2n  2n  ^  ^  gn  n2  1  ^gn"l  / .  N 

B  (out)  -  Li  -  v“  i-  f  -  1  -  ■ "  -  —77  g  +  —  (A-20 ) 

n  '  K  r  dr  2  n  >.  2  2  &n  r  or  v  / 

r  or  r  J 

.  _  rftf  2  r*g  - 

c  (out)  ,  A±2H.  n  .2  g  .  «L  f  .nJll 
n  '  r  |^or  r  cn  r  n  or  J 

rd2f  dg  ng  -j 

,  n  n  °n  °n 

+  x  br  *  7  sr  -  -s\ 


Again ;  using  the  equations  for  f^  and  g^  given  earlier  we  can  write  the  above 
coefficients  finally  as 

00 

(out)  =  i J'  F1‘'l’:?(§1)  cosh  nu1  "  +  2  cosh2  u.J  du^^ 


-  J  G^vSg)  Sinh  2Ug  sinh  nUg  du2 

C2  ° 


f  i  \n  p 

B^out)  =•  J  Fn  sinh  2u^  sinh  nu^  du1 

C^  ^  0 


'-l^  r  n-2  ■ 

J  G  'Ig.1  cosh  nu2  cosh  2ug  dug 


c-2  J0~  -2- 


(A-21) 


C(out)  =  cosh  nu  —  -  2  sinh2u  dun 

c"  J  0  1-l  1 


f  J  Gn  tig;  sinh  2u2  sinh  nUg  dUg 

C2  ° 
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APPENDIX  B 
STATIC  SOLUTION 


B.l  Dilatational  Wave 

The  static  solution  presented  here  is  based  on  the  application  of 
the  theory  of  stress  functions  presented  oy  Timoshenko  (9) .  Under  static 
conditions ,  it  can  be  assumed  that  at  large  distances  from  the  boundary  of 
the  cavity  the  state  of  stress  in  the  medium  is  equal  to  the  stress  field  in 
the  medium  behind  the  front.  In  polar  coordinates  this  is 


rr 


■ 


a  cos  20 
P 


ar& 

°9d 


=  (ttt)  ap  sin  20 

"  ^  °p  +  °p  cos  29 


(B-l) 


is  the  stress  in  the  medium  in  the  direction  of  wave  propagation  and  va 
is  the  stress  parallel  to  the  wave  front.  v>  derived  from  the  assumption 
that  there  is  no  strain  parallel  to  the  wave  front,  is  equal  to 


V  =  jJL.  (B-2) 

It  can  be  seen  from  equations  (B-l)  that  the  n  =  0  and  2  modes  describe 
exactly  the  free  field  stresses.  Therefore,  the  unknown  stresses  at  the  edge 
of  the  cavity  can  likewise  be  expressed  in  terms  of  these  modes 


a 


rr 


Pq  -  P2  cos  20 


(B-3) 


V  =  S2  8in  29 

where  PQ,  an(i  ^2  are  unicnown  mCKial  components  of  stress  acting  at 
the  boundary.  These  same  stresses  must  act  on  the  shell.  Thus  the  boundary 
stresses  in  the  medium  and  shell  can  be  illustrated  as  shown  in  Fig.  9. 
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The  equations  of  equilibrium  in  the  medium  under  plane  strain 
conditions  are  satisfied  by  the  following  expressions  for  the  components  of 
stress 


o 


rr 


i  an  i  a2n 
=  r2^2 


_  _1  dfl  1 

rG  2  §0  "  r  drcte 
r 


£a 

Sr2 


(B-4) 


where  ft  is  the  stress  function  in  terms  of  r  and  G. 


B.ll  n  =  0  Mode 

The  general  solution  of  the  stress  function  for  the  n  =  0  mode  is 
given  by  Timoshenko  as 

2 

ft  *  K  log  r  +  Mr 

from  which  can  be  derived  the  stress  components 


a 

rr 


o 


rQ 


=  0 


(B-5) 


<J 


QQ 


V  SM 

r 


The  coefficients  K  and  M  are  determined  from  the  states  of  stress  at  r  -4  « 
and  r  =  R 


(B-6) 


For  the  conditions  of  plane  strain,  the  strains  are 
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-r  I|0rr-vaeel 


■it 

it 


aQQ  -  v  arr 


(B-7) 


”  ipi  °r9 


where 


E  = 


E 


1  -  V 


modulus  of  plane  strain 
for  the  medium 


Displacement  components  are  found  by  suitable  integration  of  the  following 
equations  for  the  strains  in  the  medium 

du 

er  =  ¥ 


u  ,1  dv 
€e  =  r  r  Se 


y  =  1  8u  ,  c)v  _  V 

' tQ  r  Ss  Sr  "  r 


(B-8) 


Displacements  in  the  medium  for  the  zero  mode  are 


u 

R 


v 

R 


0 


(B-9) 


The  corresponding  displacement  components  of  the  shell  with  an  exterior 
compressive  force  PQ  are  found  in  Flugge  (2)  to  be 


where 


(B-10) 


modulus  of  plane  strain 
for  the  shell 


The  unknown  stress  Pq  is  now  determined  by  equating  the  radial  displacement 
components,  to  get 


P 

o 


(l+v)  a 


B  R  .  — " 

■r-  •  r  +  1  +  V 

J!  A  VJ 


(B-ll) 


After  Pq  is  determined,  stresses  in  the  medium  at  any  radius  can  be  found 
using  equations  (B-5).  Equation  (B-10)  gives  the  displacement  of  the  shell; 
and  the  hoop  and  bending  stresses  of  the  shell  are,  respectively 


<*09 


P  R 

o 

A 


(B-12) 


B.12  n  =  2  Mode 


The  general  solution  in  terms  of  a  stress  function  ft  is  given  as 


ft  =  |^kr^  +  -4^  +  Ifj  cos  20 


from  which  the  stress  components  in  the  medium  become 


rr 


,  .  [at  ♦  *  *  Ml  cos  „ 

L  r  rdj 
=  j^2K  +  cos  29 

-  [a  -  nj  -  “1  sin  20 
'*  r  r  ■* 


r0 


(B-13) 


Coefficients  are  determined  as  before  from  the  states  of  stress  at  r 
and  r  =  R 


K  =  (^)  °r 


(B-14) 
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m  R 
M  =  T 


2iizxl  Op  -  P2  -  2S 


J 


R2  r 

N  =  — 

1  2 


_.(l-v)  ap  +  P2  +  sj 
The  displacement  components  of  the  medium  at  the  boundary  can  be  expressed 


as 


u 

R 


v 

R 


E 


2(2-v)S, 


2(v-l)  +  P2  +  5 


1 


cos  29 


—  2(l-v)  -  P2  -  Sgj  sin  20 


(B-15) 


and  the  corresponding  components  for  the  shell  are 


.  r  %  g .  %  _e£i  cos  2e 

R  I.  E  91  e  J 

s  s 

Vs  [P2  /  R5  R  x  S2  ,  R5  R  xl  a 

t  =  Lr ( nr +  sa  }  +  r ( w +  sin  2® 


(B-16) 


The  unknown  boundary  stresses  and  are  determined  by  equating  the  dis¬ 
placement  of  the  shell  and  medium  at  the  boundary.  The  hoop  and  bending 
stresses  in  the  shell  are  then  found  from  the  following  equations 


PpH  2S  R 

°00  =  (jr +  —}  cos  29 


°SB  ■  hi  <2P-  *  s-> c°s  29 


(B-17) 


Displacements  of  the  shell  and  stresses  in  the  medium  can  be  determined  using 
equations  given  earlier  in  the  discussion.  Total  static  stresses  and  dis¬ 
placements  are  merely  the  sum  of  the  modal  values  for  any  angle. 
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B.2  Shear  Wave 

The  free  field  shear  wave  stresses  expressed  in  polar  coordinates 
include  only  the  n  =  2  mode 


a  =  a  sin  20 
rr  s 

ar0  =  as  cos  20  (B-l8) 

°00  =  -os  sin  20 

where  a  is  the  amplitude  of  the  incoming  wave*  The  static  solution  is 
s 

obtained  exactly  as  in  the  case  of  the  dilatational  wave  except  for  the 
interchange  of  sines  and  cosines  resulting  from  the  difference  in  geometry. 
The  resulting  equations  are 

a.  Displacement  components 


saB?] 

18  E  iJ 

s 


sin  20 


cos  20 


b.  Shell  Stresses 


(B-19) 


o 


o 


tp2R 

00  - 

"  |> 

sB  = 

-  i 

121 

2S_R"| 


P2  +  S2)  sin  20 


(B-20) 


c. 


Medium  Stresses 

°rr  =  "  ^  +  %  +  ”*§)  sin  20 
r  r 

or0  =  -  (2A  -  ^  cos  20 

r  r 

aQQ  =  (2A  +  sin  20 


r 


(B-21) 
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where 


a 


In  the  machine  program,  the  n  =  0  mode  is  purely  rotational  and 
thus  we  get  a  static  value  of  tangential  displacement.  This  is  calculated 
as  follows 


9  «i  7  =iJL 

vot  2  xy  2  ji 


.  a 

v  _  Q  _  1  s 
R  rot  2  n 


(B-22) 
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DILAIATIOKAL  WAVE 


tine  * 


Dev.  Fro* 

Nod*  Quant*  Stat.  Nach.  Static 


Dtr.  Froi 


0 

-d. 

A  *000 

4.006 

.15* 

0 

-Oh 

4.000 

4.006 

.15 

0 

<r*w 

.100 

.100 

.00 

Node 

Quant 

Stat. 

Nach. 

Static 

2 

-OL 

4.613 

4.599 

•  30* 

2 

fl 

3.845 

3.856 

.29 

2 

Off 

3.077 

3.114 

1.20 

2 

<r»t> 

.346 

•  345 

.29 

*  Stresses  flvea  1b  units  of  |C>I 
Displacements  defined  by  I41.  (2-90)  and  (2-36) 


SOAR  NAVI 

tine  c  20  1/ 

Node 

Quant. 

Stat.  Nach. 

Dev.  Fro* 
8tatlc 

0 

0 

5.000  4.996 

.08* 

2 

OL 

13.840  13.930 

.65 

2 

11.535  11.539 

.03 

2 

-<*• 

9.232  9*149 

.90 

2 

(Tib 

1.038  1.045 

.67 

TABIX  1 

OQNFAHSOS  OF  SBOIC  AID  NAGUIB  VALUES  FOR  NODAL  RESPONSE  OF  TKK  SELL 
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DILATATXOIAL  WAVS 


tin*  *  20  R/ci 


1.0  R 

warn 

warn 

*.0  R 

Mod* 

Quant? 

3Ut. 

Mach. 

— — - - 

SUt.  Mach. 

Stat. 

Mach. 

SUt. 

Mach. 

0 

-  OVr 

.200 

.200 

•  550 

.551 

.615 

.616 

.638 

.639 

0 

-OW 

1.133 

1.135 

.783 

.785 

.718 

.720 

.696 

.697 

2 

■  OVr 

-.153 

-.155 

.111 

.112 

.226 

.229 

.271 

.273 

2 

ov# 

.307 

.312 

.kse 

.*29 

•  38k 

.386 

.363 

.366 

2 

0»» 

.872 

A£A 
•  OOO 

.367 

.368 

•  3*0 

•  3»* 

.335 

.337 

*  8tniMi  glT*n  In  unit*  of  ld*| 


SIAR  WAVS  tiM  s  20  R/e2 


1.0  R 

2.0  R 

3.0  R 

*.0  R 

Mod* 

(feiant. 

SUt.  Mach. 

BUt. 

Mach. 

BUt. 

Mach. 

SUt.  Mach. 

2 

OW 

-.*60  -.A  5  5 

.332 

.331 

.675 

.678 

.81*  .815 

2 

0V# 

•921  .913 

1.283 

1.286 

1.151 

1.15* 

1.090  1.092 

2 

-0W 

2.616  2.636 

1 _ 

1.101 

1.10* 

1.020 

1.022 

1.006  1.009 

XABLB  2 


ooMMRiaoa  or  bbatic  ard  macurb  valubs  for  nodal  sirbssss  n  m  medium 
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\ 

«».{50') 

r  ~i 

rar* 

xj 

CTa*^ 

429!) _ 

nr 

*i 

Static 

Rax. 

Static 

Max. 

.01 

-9.21 

-10.12 

1.10 

3.12 

-2.49 

-2.73 

1.09 

3.13 

.05 

-7.08 

-  7.74 

1.09 

2.88 

-2.01 

-2.19 

1.09 

2.88 

.10 

-5.56 

-  6.15 

1.10 

2.25 

-1.67 

-1.84 

1.11 

2.25 

.15 

-4.60 

-  5.21 

1.13 

2.25 

-1.45 

-1.65 

1.13 

2.25 

.20 

-3.9^ 

-  4.52 

1.14 

2.25 

-1.30 

-1.50 

1.15 

2.25 

Tt 

<J»b  (90*) 

nr 

ta 

^rjrdSCL) 

nr 

m 

Static 

Max. 

.  .  . 

[  Static  | 

Max 

.01 

.10 

.11 

1.11 

3.38 

-  .09 

-  .10 

1.10 

3.25 

.05 

.45 

.50 

1.11 

3.25 

-  .35 

-  .39 

1.11 

2.88 

.10 

.83 

.94 

1.13 

3.00 

-  .55 

-  .63 

1.14 

2.50 

.15 

1.18 

1.36 

1.15 

2.88 

-  .68 

-  -79 

1.17 

2.50 

.20 

1.50 

1.76 

1.18 

2.88 

-  .75 

-  .91 

1.21 

2.50 

It 

OL 

Jol _ 

nr 

ta 

$ 

(45*) 

nr 

ta 

Static 

Max. 

Static 

Max. 

.01 

-9.65 

-10.64 

1.10 

3.12 

4.67 

5.23 

1.12 

3.75 

.05 

-8.61 

-  9.50 

1.10 

3.00 

3.85 

M3 

1.13 

3.50 

.10 

-7.76 

-  8.62 

1.11 

2.76 

3.33 

3.82 

1.15 

3.25 

.15 

-7.15 

-  8.01 

1.12 

2.88 

3.02 

3.55 

1.18 

3.25 

.20 

-6.66 

-  7.56 

1.13 

2.88 

2.80 

3.37 

1.20 

3.38 

*  Ratio  of  Mu.  to  Static  #  Tim  1.  transit  tlass  at  which  Max.  occurs 

TABU  3 


ORCT  or  RELATIVE  1BICK1BSS  CT  SHE 
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0m» 

(90* ) 

DU* 

M 

O«iw.(90*) 

H 

■ 

8tatlc 

Max. 

Static 

Max. 

2 

-  4.16 

-  4.57 

1.10 

2.76 

-2.27 

-2.51 

l.il 

2.76 

4 

-  7.08 

-  7.74 

1.09 

2.88 

•2.01 

-2.19 

1.09 

2.88 

6 

-  9-33 

-10.11 

1.08 

2.76 

-1.81 

-1.97 

1.08 

2.76 

10 

-12.57 

-13.41 

1.07 

2.76 

-1.55 

-1.65 

1.07 

2.76 

16 

-15.72 

-16.47 

1.05 

2.88 

-1.31 

-1.37 

1.05 

2.76 

\ 

(*>") 

DU 

ta 

(90*) 

DU 

ta 

Static 

1 

• 

Static 

Max. 

2 

.23 

.26 

1.12 

3.25 

-  .21 

-  .23 

1.13 

3.12 

4 

.45 

.50 

1.11 

3.25 

-  .35 

-  .39 

1.11 

2.88 

6 

.64 

.71 

1.11 

3.12 

-  .47 

-  .51 

1.10 

2.76 

10 

1.01 

1.12 

1.10 

3.12 

-  .63 

-  .68 

1.08 

2.62 

16 

1.53 

1.68 

1.09 

3.12 

-  .78 

•  .83 

1.06 

2.62 

*  (0*) . 

a 

(45*) 

If 

Static  1  Max. 

1 _ 

DU 

ta 

Static 

Max. 

DU 

ta 

2 

-  4.60 

-  5.12 

1.11 

3.00 

2.14 

2.43 

1.13 

3.62 

4 

-  8.61 

-  9.50 

1.10 

3.00 

3.05 

4.33 

1.13 

3.50 

6 

-12.22 

-13.37 

1.09 

2.88 

5.32 

5.96 

*.12 

3.38 

10 

-18.65 

-20.20 

1.08 

2.88 

7.92 

8.82 

1.11 

3.25 

16 

-27.11 

-29-15 

1.08 

3.00 

11.41 

12.64 

1.11 

3.12 

*Batio  of  Nix.  to  Static  ^Tlae  in  transit  tlaes  at  which  Max.  occurs 

TABLE  4 


EFFECT  OF  RATIO  CT  MODULI  GT  ELASTICITY 


> 

■E9I 

1 

nr* 

H 

■E9R5M 

OUT 

■ 

Static 

Msx. 

Static 

Msx. 

1 

-7.08 

-7.67 

1.08 

3.00 

-2.01 

-2.17 

1.08 

2.88 

3 

ft 

-7.7k 

1.09 

2.88 

If 

-2.19 

1.09 

2.88 

5 

ft 

-7.97 

1.13 

2.50 

ft 

-2.26 

1.12 

2.50 

7 

ft 

-8.30 

1.17 

2.61 

♦f 

-2.3*i 

1.17 

2.62 

10 

ft 

-8.70 

1.23 

2.88 

»» 

•2.45 

1.22 

2.88 

V 

<r»t 

wsm 

H 

■ 

mm 

nr 

- 1 

i 

ta 

Static 

Msx. 

Static 

Max. 

1 

.45 

.49 

1.10 

3.25 

-  .35 

-  .33 

1.09 

3.00 

3 

N 

.30 

l.u 

3.25 

ft 

-  .39 

l.U 

2.88 

5 

m 

.50 

1.13 

3.12 

N 

-  .*>1 

1.15 

2.88 

7 

t« 

.51 

1.15 

3.25 

n 

-  .43 

1.21 

3.12 

10 

M 

.53 

1.18 

3.50 

*t 

-  .46 

1.30 

3.50 

V 

OL 

JSL) _ 

nr 

ta 

!  * 

(45*1 

nr 

*i 

Static 

|  Static 

Max. 

1 

-8.61 

-9.38 

1.09 

3.00 

3.85 

4.27 

l.U 

3.50 

3 

ft 

-9.50 

1.10 

3.00 

ft 

4.33 

1.13 

3.50 

5 

ft 

-9.64 

1.12 

2.88 

ft 

4.41 

1.15 

3.62 

7 

ft 

-9*77 

i.i4 

3.00 

ft 

4.45 

1.16 

2.76 

10 

ft 

-9.95 

1.16 

2.62 

if 

4.72 

1.23 

3.00 

*ktto  of  Hue.  to  Static  la  transit  time  at  which  Max.  occurs 

TABLE  5 


VPICT  or  MASS  RATIO 


Static 


DU*  |  t«# 

[  3tatlc  [ 

Max.  | 

..10 

3.25 

•2.08 

-2.29  1 

-°9 

3.12 

•2.06 

-2.2k  1 

..09 

3.00 

•2.02 

-2.20  1 

l.09 

2.88 

•2.01 

-2.19  1 

.17 

3.12 

-1.91 

•2.22  1 

Otb 

Static 

Max. 

oc 

±0*1.  _ 

Static 

Max. 

Static 


3HH 


0 

-10.33 

-11.48 

1.11 

3.25 

6.00 

6.9* 

1.16 

3.38 

.1 

-  9.88 

-10.90 

1.10 

3.12 

5.31 

6.08 

1.14 

3.38 

.2 

-  9.12 

-10.03 

1.10 

3.00 

4.4o 

4.98 

1.13 

3.38 

.25 

-  8.61 

-  9.50 

1.10 

3.00 

3.85 

*.33 

1.13 

3.50 

.k 

-  6.54 

-  7.61 

1.16 

3.50 

1.77 

1.97 

1.11 

*.38 

*Ratlo  of  Max.  to  Static 


^Ti a*  In  transit  tlaes  at  which  Max.  occurs 
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1.2 
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PIO.  26  RADIAL  STRESS  IE  THE  MEDIUM  AT  THE  BOUHDAR  FOR  VARIOUS  THU 


1.2 
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FIO.  31  COMPARISON  OF  SHORT-TD®  RESULTS  FOR  IECUEET  P  WAVE 
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FIG.  33  COMPARISON  OF  SHORT* TIME  RESULTS  FOR  HICIEEHT  P  WAVE 
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H 
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i r\ 
i 
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no 

i 


P  WAVE 
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